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Kurzfassung 



In dieser Arbeit prasentieren wir das Computeralgebra Paket HarmonicSimis und 
dessen theoretischen Hintergrund fiir die Verarbeitung von harmonischen Summen und 
damit verbundenen Objekten wie zum Beispiel Euler-Zagier Summen und harmonis- 
che Polylogarthmcn. Harmonische Summen und vcrallgemeinerte harmonische Sum- 
men treten als Spezialfalle der sogenannten d'Alembert Losungen von Rekurrenzen 
auf. Wir zeigen, dass die harmonischen Summen eine quasi-shuffle Algebra bilden und 
prasentieren eine Methodc zum Auffinden algebraisch unabhangiger harmonischer Sum- 
men. Aui3erdem definieren wir eine Differentiation auf den harmonischen Summen iiber 
eine Erweiterung der Melhn Transformation. Dabei treten neue Relationen zwischen 
den harmonischen Summen auf. Zusatzlich stellen wir einen Algorithmus vor, der be- 
stimmte verschachtelte Summen in Ausdriicke bestehend aus harmonischen Summen 
umschreibt. Durch nicht triviale Beispiele illustrieren wir, wie diese Algorithmen in 
Zusammenarbeit mit dem Summationspaket Sigma die Berechnung von Feynmaninte- 
gralen unterstiitzen. 
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Abstract 



In this work we present the computer algebra package HarmonicSums and its theoretical 
background for the manipulation of harmonic sums and some related quantities as 
for example Euler-Zagier sums and harmonic polylogarithms. Harmonic sums and 
generalized harmonic sums emerge as special cases of so-called d'Alembertian solutions 
of recurrence relations. We show that harmonic sums form a quasi-shuffle algebra 
and describe a method how we can find algebraically independent harmonic sums. In 
addition, we define a differentiation on harmonic sums via an extended version of the 
Mellin transform. Along with that, new relations between harmonic sums will arise. 
Furthermore, we present an algorithm which rewrites certain types of nested sums into 
expressions in terms of harmonic sums. We illustrate by nontrivial examples how these 
algorithms in cooperation with the summation package Sigma support the evaluation 
of Feynman integrals. 
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Notations 



N N = {1, 2, . . .}, natural numbers 

Z integers 

Q rational numbers 

R real numbers 

Sai,a2,...i''^) harmonic sum; see page 3 

Zai,a2,-:i^) Euler-Zagier sum; see page 5 

Hmi,m2,...(^) harmonic polylogarithm; see page 29 

Sp{n) the set of polynomials in the harmonic sums with integer indices; see page 

S{n) the set of polynomials in the harmonic sums; see page 15 

Ip an ideal of Sp{n); see page 14 

I an ideal of 5(n); see page 15 

Lu the shuffle product; see pages 8 and 18 

* the quasi-shuffle product; see page 9 

A* the free monoid over A; see page 7 

A'^ the set of non-empty words over A; see page 8 

R {A) the set of non-commutative polynomials over R; see page 7 

Lyndon(A) the Lyndon words over A; see page 8 

£() i{ w) gives the length of a word w 

In the number of Lyndon words; see page 16 

sign() sign(a) gives the sign of the number a 

II \w\ gives the degree of a word e; |a| is the absolute value of the number a 
A a A 6 = sign(a)sign(5)(|a| + |6|); see page 4 

/x(n) the Mobius function; see 16 

(^n Cn is the value of the zeta-function at n 

part (a) see page 61 

rev(a) see page 61 

1„ a vector of length n where all components are 1 

On a vector of length n where all components are 

M {f{x), n) the Mellin transform of f(x); see page 43 

M"*" (f(x),n) the extended Mellin transform of f(x); see page 43 
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Chapter 1 



Introduction 



Multiple harmonic sums, see e.g., [5, 41] and Nielsen-type integrals, associated to them 
by a Mellin transform, emerge in perturbative calculations of massless or massive single 
scale problems in quantum field theory; for more details how these sums emerge we refer, 
i.e., to [8]. Due to the complexity of higher order calculations the knowledge of as many 
as possible relations between the finite harmonic sums is of importance to simplify the 
calculations and to obtain as compact as possible analytic results. The multiple finite 
harmonic sums 5'ai,a2,...,afc('^) defined as 

Sa. aM= E ^^"^f ■ ■ ■ "'''"^f (1-1) 

n>ii>i2>--->?fc>l *1 'fc 

yield a very convenient description for these field theoretic quantities [6]. 
Harmonic sums form a quasi-shuffie algebra [19, 20, 21] and the algebraic relations of 
the harmonic sums are well known [5]. The number of the respective basis elements 
are counted by the Lyndon words. All further relations between the nested harmonic 
sums are called structural relations. They result from the mathematical structure of 
these objects beyond that given by their indices [6, 8, 7]. 

In physical applications [11, 4, 5, 12, 6] an analytic continuation of the nested harmonic 
sums is required and eventually one has to derive the complex analysis for nested 
harmonic sums. In passing, n takes values n G Q and n G M, which leads to new 
relations [6]. Considering n G Q leads, e.g., to half-integer relations, while considering 
n G M will allow us, e.g., to differentiate with respect to n and so new structural 
relations will arise. 

Multiple harmonic sums can be represented in terms of Mellin integrals [41, 11], in fact 
they are the Mellin transforms of harmonic polylogarithms [33], which belong to the 
Poincare iterated integrals see [31, 24]. This fact allows us to establish a differentiation 
on the harmonic sums. 

The package Harmonic Sums tries to combine all these features and generalizes some 
aspects to Euler-Zagier sums and S-sums (see [28]). Besides basic operations such 
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as multiplication and evaluation of harmonic sums, harmonic polylogarithms, Eulcr- 
Zagicr sums and S-sums it provides procedures to compute the Mellin and the inverse 
Mellin transform, to differentiate harmonic sums and to find algebraic and structural 
relations. Tables of these relations are provided and procedures to apply these relations 
are included in the package HarmonicSvims . 

With the package Sigma [36] we can simplify nested sums such that the nested depth 
is optimal and the degree of the denominator is minimal. This is possible due to a 
refined summation theory [39, 37] of HS-fields [22]. The package HarmonicSums finds 
sum representations of such simplified nested sums in terms of harmonic sums as much 
as it is possible. 

The remainder of the thesis is structured as follows: In Chapter 2 harmonic sums and 
similar strTictures are introduced. We show that harmonic sums form a quasi-shuffle 
algebra and we show the connection to Lyndon words. Finally, we outline a way how 
we can find algebraic relations between harmonic sums and how wc can apply these 
relations. While Chapter 2 only deals with algebraic relations we introduce harmonic 
polylogarithms in Chapter 3, which leads to differentiation of harmonic sums via 
Mellin-transform. Along with that new relations, structural relations will pop up. 
Chapter 4 deals with half-integer relations between harmonic sums and in Chapter 5 
we show how a certain kind of nested sums can be rewritten in terms of harmonic sums. 
In the last chapter we give a nontrivial example to illustrate how these algorithms in 
cooperation with the summation package Sigma support the evaluation of Feynman 
integrals. 

Acknowledgments. This diploma thesis was supported by the SFB grant F1305 and 
the grant P20347-N18 of the Austrian FWF. 

I would like to thank C. Schneider for supervising this work and his support. In 
addition I would like to thank J. Bliimlein for his helpful suggestions and comments. 
In particular it was an exciting experience to be invited at DESY, Zeuthen. 
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Chapter 2 

Algebraic Relations between 
Multiple Harmonic Sums 

In this section we will define harmonic sums and relate to them the Euler-Zagier sums. 
In the following we will rely on the fact that harmonic sums form a quasi-shufRe algebra. 
Quasi-shuffle algebras are connected to the free polynomial algebra on the Lyndon 
words. It will turn out that the number of Lyndon words allows one to determine 
the number of algebraically independent harmonic sums (in the context of quasi-shufSe 
algebras) and gives at least an upper bound for the number of algebraically independent 
harmonic sums (considered as sequences). In the end of this chapter we will present 
a way how we can derive these algebraically independent harmonic sums following the 
ideas from [5]. 



2.1 Definition and Product of Harmonic Sums 



We start by defining multiple harmonic sums; see, e.g., [5, 41]. 
Definition 2.1.1. ForkeN,neZ and Ui G Z/{0} with 1 < i < k, 



S{n) 



1, n > 
0, n < 0, 

sign(ai) 



1=1 



k is called the depth and = |ai| H h |afe| is called the weight of the harmonic sum 

Sai,...,ak{iT')- Note that an equivalent definition is given in (1.1). 

Notation 2.1.2. Sometimes we will write Sa-i^a^as-'-aki^) instead of Sai,a2,a3,...,ak{''^)- 
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+ 



Figure 2.1: Sketch of the proof of the multiphcation of harmonic sums (compare [28]). 



As we will see later in Section 2.3 harmonic sums obey an algebra. In this algebra the 
following property is crucial: a product of two harmonic sums with the same upper 
summation limit can be written in terms of single harmonic sums: for n € N, 

Sau...,ak (n) * Sb^,...,bt {n) = ^ ^^^^17^ '5'a2,...,afc (0 * Sbi,...fy (?) 

i=l ^ 



i=l 
n 



sign(Qi * biY 
2^ ,|ai|+|bi| S<^2,-MV * Sb,,...MW (2-1) 



ikil+lfeil 
1=1 

here * means the usual multiplication. The proof of equation (2.1) follows immediately 
from 

n n n i n j n 

= + " ^^^^ 

i=l j=l 1=1 j=l J=l '=1 i=l 

for a graphical illustration of the proof see Figure 2.1. As a consequence, any product 
of harmonic sums can be written as a linear combination of harmonic sums by iterative 
application of (2.1). We give an example for the product of two harmonic sums. 

Example 2.1.3. For n>0, 

Saiin) * Sbi,...,biin) = S'ai,6i,...,6,(n) + Sbi,ai,b2,-Mi'^) ^ Sb^fi2,-,bi,ai{n) 

-'S'aiA6i,fe2,...,6(H Sb^,b2,...,aiAbm(.'ri)- (2-2) 



Here the symbol A is defined as 

a Ab = sign(a)sign(6)(|a| + |6|). 



Using the package HarmonicSums in Mathematica we can do this automatically: 

HarmonicSums session. 
in[i]— S[l,4,n] S[2,-3,n]//LinearExpand 
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Out[i]= S[3,-7,n] - S[l,2,-7,n] - S[l,6,-3,n] - S[2,-4,4,n] - S[2,l,-7,n] - S[3,-3,4,n] - 
S[3, 4, -3, n]+S[l, 2, -3, 4, n]+S[l, 2, 4, -3, n]+S[l, 4, 2, -3, n]+S[2, -3, 1, 4, n]+S[2, 1, -3, 4, n] + 
S[2,l,4,-3,nl 



Notation 2.1.4. iuras out later, the following modified notation (also used in 

[4IJ) for multiple harmonic sums will be useful. In this notation we will only allow the 
indices 1, 0, and—1. An index zero which is followed by an index that is nonzero means 
that 1 should be added to the absolute value of the nonzero index. So for example: 

S'0,0,1,0 -1,0,0,0,1 H = 'S'3 ,-2,4(n) 
We will refer to this modified notation as the expanded version. 



Using the modified notation it is quite easy to calculate the number N{w) of multiple 
harmonic sums for a given weight w. Since the length of the index set in the expanded 
version is exactly the weight of the harmonic sum, we have to consider all possible index 
sets of length w. Each index except the last can take the three values 1, 0, and — 1. 
The last index can just be 1 or -1. Hence we get the following formula, compare [41] 
or [5]: 

iVH = 2-3"'-^ 



2.2 Euler-Zagier Sums 



Now we define in a similar way the so called Euler-Zagier sums, see [45, 17] or compare 

[281: 



Definition 2.2.1. For k eN, n eZ and ai € Z/{0} with l<i<k, 

Z{n) = 



1, > 

0, n < 0, 



Za,,...,akin) = 2^ Z„,,...,„^(z-1). 

i=l ^ 



Again k is called the depth and id = |ai| + • • • + jajtl is called the weight of the Euler- 
Za gier sum ^ai,...,afc('^)- An equivalent definition is given by 

Za,,...,aAn) - 2^ . [2.6) 

n>ii>i2>-->ik>l ^1 ^fe 
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The product for these sums can be expanded with the formula 

^ai,...,afc(n)*^6i,...,6i(n) = Yl ^'^^M -^°2.-,afc (0 * ^bi,...,bi (^) 

i=l ^ 

Asign(6i)^ 

i=l ^ 

A sign(ai * 61)* 
1=1 

Again we can perform this expansion using the package Harmonic Sums: 



HarmonicSums session. 






in[2]:= Z[l,4,n] Z[2,-3,n] 






Out[2]= Z[3,-7,n] + Z[l,2, -7,n] + Z[l,( 


3,-3,7-1] + Z[2,-4,4,n] + Z[2,l,- 


7,n] + Z[3, -3,4, ri] + 


Z[3,4,-3,n] + Z[l,2,-3,4,n] + 


Z[l,2,4,-3,n] + Z[l,4,2,-3,n] 


+ Z[2,-3, l,4,n] + 


Z[2,l,-3,4,n] + Z[2,l,4,-3,n] 







Multiple harmonic sums and Euler-Zagier sums are closely related. In particular, we 
can easily convert them from one representation to the other one, compare [28]. Namely 
by using 

^ sign(ai)^i sign(a2)'^ 

il=l ^1 t2=l *2 

+S'aiAa2,a3,--('^)! 

we can recursively convert multiple harmonic sums to Euler-Zagier sums. Similarly one 
can exploit 

^ sign(ai)^i ^ sign(a2)^^ ^ 

Za„...{n) = X. .\a,\ 2^ .\a,\ ^aa,..., (^2 - 1) 

41 = 1 ^1 12=1 *2 

~-^aiAa2,a3,...(?^)) 

in order to convert Euler-Zagier sums to multiple harmonic sums. 
Example 2.2.2. For n G N, 



Si,3A{n) = Mn) + Zijin) + Z4,4{n) + ^i,3,4(n), 
Zi,3,4{n) = Ssin) - Sij{n) - S4,4{n) + S'i,3,4(n). 
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The package HarmonicSums provides procedures for these conversions: 
HarmonicSums session. 

In[3]:= ZToS[Z[l, 3, 4, n]] 

Out[3i= S[8, n] - S[l, 7, n] - S[4, 4, n] + S[l, 3, 4, n] 
ln[4]:= SToZ[S[l, 3, 4, n]] 

Out[4]= Z[8, n] + Z[l, 7, n] + Z[4, 4, n] + Z[l, 3, 4, n] 



2.3 The Shuffle and the Quasi-shuffle Algebra 



In Section 2.4 we present the algebra of multiple harmonic sums, but first we have to 
give some basic definitions; in particular, we will introduce the shuffle and the quasi- 
shuffle algebra. 

Definition 2.3.1 (Graded algebra). A graded algebra over a ring-' R is an R-algebra 
E that has a direct sum decomposition 

E = ^Ei = Eo®Ei®E2®--- 

i 

such that for any x G Ei, and y G Ej we have xy G -Ej+j, i.e., EiEj C Ei^j and such 
that for any r & R and x & Ei we have rx G Ei, i.e., REi C Ei. 

Definition 2.3.2 (Non-commutative Polynomial Algebra). Let A be a totally ordered, 
graded set. The degree of a G A is denoted by \a\ . Let A* denote the free monoid over 
A, i.e., 

A* = {ai-- ■ an\ai e A,n>l}\J {e} . 

We extend the degree function to A* by \a1a2 • • • an| = |ai| + |fl2| + • • • + |on| for ai G A 
and \e\ = 0. Let RD Q be a commutative ring. The set of non- commutative polynomials 
over R is defined as 



R{A) := J '^^'^ 



ry, E R,ry, = for almost all w 



^Throughout this thesis we assume that all rings and fields contain Q. The ground field, denoted by 
K and the ground ring, denoted by R, are always commutative. 
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Addition in R (A) is defined component wise and multiplication is defined by 

^ a-wW ■ ^ byjW := ^( ^ aabv)w. 

■WW w uv=w 

Remark 2.3.3. We refer to elements of A as letters and to elements of A* as words. 
We call e the empty word. The length of a word w = 0102 . . . is denoted by i{w) = k. 
The set of non-empty words is denoted by A"^. 

Definition 2.3.4 (Shuffle algebra). On A*, the free monoid over A, we define the 
shuffle product as follows. If 

w = oi 02 . . . Ofc and v = bib2 ■ ■ - bi 

then 

uiue = u, 

uiuv = ai ■ (uiiuv) + bi ■ (uiuvi). 
This product is extended to R{A) by linearity. 

Remark 2.3.5. R is an easy exercise to verify that {R{A) , lu) is an associative, com- 
mutative and graded K-algebra, which is in the following called the shuffle algebra. 

.1 Lydon words 

Suppose now that the set A of letters is totally ordered by <. We extend this order to 
words with the following lexicographic order: 

M < w for u G and A* , 

WiaiW-2 < Wi 02103, ifai < 02 for 01,02 G A and wi,W2,W'i € A*. 

Definition 2.3.6. For a word w = pxs with p,x, s E A* , p is called a prefix of w if 
p ^ e and any s is called a suffix of w if s ^ e. A word w is called a Lyndon word if 
it is smaller than any of its suffixes. The set of all Lyndon words on A is denoted by 
Lyndon(A). 

Note that any word w ^ A* can be factorized uniquely in a lexicographically decreasing 
product of Lyndon words [3, 25, 14]. Therefore it can be written as w = ti ■ ■ ■ , 
where Ij G Lyndon(^) and > if i < j. The following theorem was first obtained by 
Radford [32] a proof can be found in [34] : 
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Theorem 2.3.7. The shuffle algebra {R{A) , lu) is the free polynomial algebra on the 
Lyndon words, i.e., {R{A) , lu) ^ (i?[Lyndon(^)], lu). 

In other words, any non-commutative polynomial in R{A), can be expressed uniquely 
as a commutative polynomial on the Lyndon words. 

Notation 2.3.8. For w ^ A* and i & N we use the abbreviation 

= wiuwiu ... luw . 

i times 

We can use the following lemma to rewrite each word as a commutative polynomial of 
Lyndon words, see [32] or [3]. 

Lemma 2.3.9. Let w = ■ ■ ■ tj^ be a word, where li, . . . ,1^ € Lyndon(A) with li > Ij 
if i < j and ii, . . . ,ik G N. The polynomial Qw = • • • ^feO'i^'^ lu • • • lu t^^'' can he 

written as w+Ry,, where Ry, contains only words whose weights are equal to the weight 
of w and which are smaller than w (with respect to our lexicographic order). 

Example 2.3.10. Let w = babab, where a < b. Obviously, w is not a Lyndon word. 
Using the lemma above recursively, we get 

w = ^binab^^ + 12aabbb + Aababb — 2abbiuab — 2biuaabb. 

2.3.2 The Quasi-Shuffle Algebra 

We define a new multiplication * on R{A) which is a generalisation of the shuffle 
product, by requiring that * distributes with the addition. We will see that this product 
can be used to describe properties of harmonic sums. In [21] (see also [19, 20]) there is a 
similar construction which for example gives rise to properties of Euler-Zagier sums. 

Definition 2.3.11 (Quasi-shuffle product). * : R{A) x R{A) — > R{A) is called 
Quasi-shuffle product, if it distributes with the addition and 

w = w * e = w, for all w G ^* , 
au*bv = a{u* v) + b{u* v) — [a,b]{u* v), for all a,b e A;u,v e A* , (2.5) 

where [■,■]: A x A A, (A = A(J {0}) is a function satisfying 

50. [a, 0] = for aU aeA; 

51. [a, b] = [b, a] for all a,beA; 

52. [[a, b],c] = [a, [b, c]] for all a, 6, c G A; and 

53. Either [a, 6] = or \[a, b]\ = \a\ + \b\ for all a,beA. 
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Note that the quasi-shuffle product gives rise to various appUcations by defining ex- 
phcit functions [•, •] accordingly. In Section 2.4 we exploit this construction to handle 
harmonic sums properties. 

Theorem 2.3.12. (compare [21]) {R{A) ,*) is a commutative graded algebra. 



Proof. It is enough to show that the operation * is commutative, associative and adds 
degrees. For commutativity we have to show wi * W2 = W2 * wi for any words Wi 
and W2. We proceed by induction on i{wi) + (-{w2). Since there is nothing to prove if 
either wi or W2 is empty, we can assume that there are letters a, b so that Wi = au and 
W2 = bv. Then (2.5) together with the induction hypothesis gives 

Wi * W2 — W2 * Wi = — [a, b]{u* v) + [b, a]{v * u), 

and the right-hand side is zero by the induction hypothesis and (SI). For associativity 
we have to prove {wi * W2) * W3 = Wi * {w2 * W3) for any words Wi,W2 and 103. 
We proceed by induction on £{wi) + i{w2) + ^(wa). Since there is nothing to prove if 
i{wi) + i{tV2) + ^(^3) < 3 or either wi, W2 or W3 is empty, we can assume that there 
are letters a,b,c so that wi = avi, W2 = bv2 and 1^3 = bv^. Then the following holds 
by using (S2), distributivity, the induction hypothesis and (2.5): 

{wi * W2) * ws = {a{vi * bv2) + b{avi * V2) — [a, b]{vi * ^2)) * CV3 

= a{{vi * bv2) * cvg) + c{a{vi * bv2) * ws) — [a, c]((wi * bv2) * V3) 
-|-6((awi * V2) * CV3) + c{b(avi * V2) * v-^) — [b, c]{{avi * V2) * W3) 

— [a, * t;2) * CV3) — c{[a,b]{vi * V2) * V3) + [[a, b], c]{{vi * V2) * vs) 
= a{vi * {bv2 * CV3)) + c{{avi * bv2) * fs) — [a, c](t;i * {bv2 * Vi)) 

+b{avi * {v2 * CV3)) — [b, c]{avi * {v2 * V3)) — [a, b]{vi * {v2 * cvs)) 
+[a, [6, c]](i;i * {v2 * V3)) 
= a{vi * {b{v2 * CV3))) + a{vi * c(6t;2 * fs)) — a{vi * [b, c]{v2 * V3)) 
+c{avi * {br}2 * V3)) — [a, c]{vi * {bv2 * V3)) + b{avi * {1)2 * cua)) 

— [b,c]{avi * {v2 * V3)) — [a,b]{vi * {v2 * CV3)) + [a, [b,c]]{vi * {v2 * V3)) 
= avi * {b{v2 * CV3) + c{bv2 * V3) — [6, c](v2 * V3)) = w\* {1112 * ws). 

To show that * adds degrees it suffices to show | wi * W2 1 = | wi | + | W2 1 • We apply 
again induction on l{wi) + £{w2). Since there is nothing to prove if either wi or W2 is 
empty, we can assume there are letters a, b so that Wi = au and W2 = bv. Recall that 
\au\ = |a| + |w| , see Definition 2.3.2. By (2.5) we get 

\wi * W2\ = \a{u* bv) + b{au* v) — [a,b]{u* v)\ . 

By (S3) and the induction hypothesis the right-hand side is smaller or equal to |wi| + 

I W2 1 . □ 
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Similar to [21] wc proceed by defining an isomorphism (p{w) : R{A) — t- R{A), see 
Theorem 2.3.15. In order to accomplish this task, we will use the following notation 
which extends the operation [•, •] as follows. Define inductively [S] G A for any finite 
sequence S of elements of A by setting [a] = a for a & A and setting [a, S] = [a, [S]] for 
any a e A and sequence S of elements of A. We get the following proposition [21]: 

Proposition 2.3.13. (i) If [S] = 0, [S'] = whenever S is a subsequence of S' . 
(ii) [S\ does not depend on the order of the elements of S. 

(Hi) For any sequences Si and S2, [Si U S2] = [[Si], [S2]] where Si U S2 denotes the 
concatenation of sequences Si and S2- 

(iv) If[S] 7^ 0, then the degree of S is the sum of the degrees of the elements of S. 

Proof, (i), (ii), (Hi), (iv) follow from (SO), (SI), (S2), (S3) respectively. □ 

A composition of a positive integer n is a sequence I = {11,12, ■ ■ ■ ,ik) oi positive integers 
such that ii + 12 + ■ ■ ■ + ik = n. We call k = the length of / and write C(n) for 
the set of compositions of n. As in [21], compositions act on words via [•, •] as follows. 
For any word w = 0102 • • • dn and composition I = {ii,i2, ■ . ■ , ik) £ C(n), set 

I[w\ = [ai, . . . ,aij[aii+i, . . . ,011+12] • • • [ai^+...+i^_^+i, . . .,an\. 

Now let (p{w) : R {A) R {A) be the linear map with (fi{e) = e and 

^{w)= Yl ■ , ■ , iii---ik)[w] (2.6) 

{h,...,ik)ec{i{w)) 

for any nonempty word w. There is an inverse i/j oi ip given by 

^{w)= J2 ^^{h...ik)[w] (2.7) 
(ii,...A)ecWi«)) ^1 

for any word w, and extended to R (A) by linearity; this follows by taking f{f) = = 
1 - = 1 - E^=o ^ = E~=i ^ in the following Lemma. {f-\t) = log(i^) = 

Eoo 
n=l n ' 

Lemma 2.3.14. (see [21]) Let f{t) = ait + a2t^ + a^t^ . . . be a function analytic at the 
origin, with ai ^ and Oj G M for all i, and let f~^{t) = bit + 62*^ + 63 . . . be the 
inverse of f. Then the map : R{A) — >■ R{A) given by 



= m ai.Oi^ ■ ■ ■ aij[w] 

i={h,.-,ik)eC{e{w)) 



(2.8) 
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for words w, and extended by linearity, has the inverse ^'^ = ^'/-i given by 

^f-i{w)= Yl hi,bi,---bij[w]. (2.9) 

Theorem 2.3.15. ip is an isomorphism of {R{A) , lu) onto {R{A) ,*) (as graded k- 
algebras). 



Proof. We follow the proof of [21, Thm. 2.5]. From the Lemma 2.3.14 </? is invertible. 
Prom Proposition 2.3.13fwJ and (2.6) it follows that |v5(tt))| < \w\. Similar, for the 
inverse we have |V'(*y)| < |^| and thus |(/:'(w)| = \w\ = |V'('"^)|- To show that is a 
homomorphism it suffices to verify that ip{vin w) = (p{v) * (^(w) for any words w, v. Let 
v= ai - ■ - am and w = hi ■ ■ ■ bn- Evidently both, ip{vuj w) and ip{v) * ip{w) are sums of 
rational multiples of terms generated by 

[SiUTi][S2UT2]---[SiUTi] (2.10) 

where the Si and Tj are sequences of ai, . . . , Om and 5i, . . . , 6„, respectively, such that 

1. for each i, at most one of the Si and Tj is empty; and 

2. the concatenation S1S2 ■ ■ ■ Si is the sequence ai,...,am and T1T2...T1 is the 
sequence 61, . . . , 6„. 

Now the term generated by (2.10) arises in (p{v) * ip{w) in only one way. The absolute 
value of the coefficient of the term given by (2.10) in ^p{v) * ^p{w) is 

1 

(length 5i)! (length ^2)! . . . (length 5/)! (length Ti)! (length Ts)! . . . (length Ti)\' 

To determine the sign of the coefficient, let V\ be the word in ip{v) and let W\ be the 
word in ip{w) such that (2.10) is a word in vi * W\. Let x, y be the lengths of W\ and 
V\ respectively. The sign is given by the following ingredients: 

• the sign of Vi in ip{v) which is (—1)"*"^; 

• the sign of W\ in <^(w) which is (— l)""^; 



the sign of (2.10) in v\ * Wi which is {-\)i^+y)-K 
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Hence the coefficient of the term produced by all possibilities of (2.10) in ip{v) * ^{w) 
is: 

(length (length Sa)! . . . (length 5/)! (length Ti)! (length Ts)! . . . (length T;)! ' 

On the other hand, the generation of a term as in (2.10) arises in ip{vujw) from 

length U Ti\ /length 5*2 U /length Si U Ti" 

length 5i J \ length 5*2 / V length Si 
(length 5i U Ti)! • • • (length 5; U T;)! 
(length Si)! . . . (length 5^! (length Ti)! . . . (length T,)! 

distinct terms of the shuffle product viuw, and after application of ip each such term 
acquires a coefficient of 

( 1 \m+n—l 



(length 5i U Ti)! • • • (length 5; U T;)! ' 
This proves that ip{viuw) = ip{v) * ip{w). □ 



It follows from Theorems 2.3.7 and 2.3.15 that {R (A) , *) is the free polynomial algebra 
on the elements {(p{w)\w is a Lyndon word}. In fact the following is true: 

Theorem 2.3.16. The quasi- shuffle algebra (i? (A) , *) is the free polynomial algebra 
on the Lyndon words, i.e., {R{A) , *) ~ (i?[Lyndon(A)], lu). 



2.4 Harmonic Sums as Quasi-Shuffle Algebra and Unique 
Representations 



Subsequently, we specialize the Quasi-shuffle algebra from Definition 2.3.11 in order to 
model the harmonic sums accordingly. We consider the alphabet 

= {1,2,3,4,...}, 

and define the degree of a letter n by |n| = n. We order the letters by n -< m for all 
n, m G N, with the usual order n < m and extend this order lexicographically. Now 
we define [a,b] = a + b for all a,b £ A and [a, 0] = for all a G Ap. This function 
obviously fulfills (S0)-(S3) and therefore {R{Ap) , *) with (2.5) is the free polynomial 
algebra on the Lyndon words by Theorem 2.3.16. The harmonic sums defined in (2.1), 
where we only consider positive integers in the index set, are completely specified by a 



Algebraic Relations between Multiple Harmonic Sums 



14 



word w & A* and the upper summation limit. E.g., the formula for the multiplication 
of harmonic sums given in (2.1) is identical to * with [a, b] = a + b: 

(n) - 5( 

a+b)(v*w) in). 

Here we need the following definition. 
Definition 2.4.1. Let a, be A*. If 

a*h = cidi H h Cmdm 

for <Zj G ^4* , Ci G M then we define 

Sa*b{n) := ciSdiX H h CmSd^in). 

Summarizing, the harmonic sums with indices in the natural numbers form a quasi 
shuffle algebra. 

We are now considering the following set: 

Sp{n) = {q{Sa„...,Sa,)\ forallrGN; ai, . . . , a, G g G M[xi, . . . , x,]} . (2.11) 

Sp{n) forms a commutative ring with infinitely many variables. In addition, we define 
the ideal Ip on Sp{n) : 

2p := {Sav{n)Sbwin) - 5'„(^6^)(ra) + Sh<^a^y,){n) - S(^aAb){v*w){n) \a,b e N; v, w € N* } . 

(2.12) 

Note that 

Sp{n)/Ip^R{A) 

by construction. In particular, the linear expansion of a G Sp{n) into harmonic sums 
gives a unique representation of a+Ip. We remark that it has been shown in [27] that 

R{A)^Sp{n)/ r-. 

where a{n) ~ b{n) ^^^k e Na{k) = b{k),i.e., the quasi-shufHe algebra is equivalent to 
the harmonic sums considered as sequences. 

We now consider the alphabet 

^ = {-1,1,-2,2,-3,3,-4,4,...}. 

We define the degree of a letter |zbn| = n for all n G N, so in each degree there are two 
letters. We order the letters by —n -< n and n ^ m for all n, m G N with n < m, and 
extend this order lexicographically. We define [a,b] = sign(a)sign(6)(|a| + \b\) = a A 6 
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for all a,b E A and [a, 0] = for all a E A. This function obviously fulfills (S0)-(S3) 
and therefore (R (A) , *) with (2.5) is the free polynomial algebra on the Lyndon words. 
The harmonic sums defined in (2.1) (now we allow negative indices) are completely 
specified by a word w e A* and the upper summation limit. 

E.g., we observe that the formula for the multiplication of harmonic sums given in (2.1) 
is identical to * with [a, b] = a Ab: 

Sav{n) * Sf,y,{n) = >S'a(„*5^) (n) + St,^av*w){''^) ~ S{aAb){v*w){''^)'i 

subsequently, we will write Sa{n)Sb{n) instead ofS'a(n) * Sb{n) for any index sets a, b. 
Summarizing, the harmonic sums with indices in Z/{0} form a quasi shufHe algebra. 

We are now considering the following set: 

S{n) = {q{Sa„...,S^)\ forallrGN; ai, . . . , a, G A*; q eR[xi, . . . ,Xr]} . (2.13) 

Moreover, we define the ideal I of 5(n) by 

(n) + Sb(av*w) (n) - S(^aAb)iv*w) (^^) \a,beZ;v,we {Z/{0})*} . 

(2.14) 

Note that 

S{n)/I^R{A) 

by construction. In particular, the linear expansion of a G S{n) into harmonic sums 
gives a unique representation of a + X. To our knowledge it has not been shown so far 
that 

R{A) ^S{n)/ ^ 

where a(n) ~ 6(n) 4^yk e Na{k) = b{k). Nevertheless, we strongly believe in this fact; 
see also Remark 2.5.1. 

Example 2.4.2. Let A = S'i,2(n)S'2(n) and B = S2infSi{n) + S2{n)S3{n) - 
S2{n)S2,i{n). Ln S{n)/X we can test equality as follows. Using the product formula 
we get 

A = -5i,4(n) - 53,2(n) + 25i,2,2(n) + -S2,i,2(n) 
B = -51,4(77,) - 5*3,2(71) + 25i,2,2(n) + S'2,i,2(n), 

thus A is equal to B. 

If we expand all products, the number of sums which are contained in the expression will 
grow relatively fast. In the previous example the expanded version of the expressions 
contains 4 sums, while there are just 2 in A. This tends to be even worse if we deal 
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with sums of higher weight or depth. Subsequently, we wiU introduee another way to 
handle the problem of unique representations such that the canonical representation of 
the expression is as small as possible (small with respect to a certain order, which has 
to be defined). We will attack this problem in Section 2.6. 



2.5 The Number of Lyndon Words or Basic Sums 



Subsequently, we will count the algebraic independent sums in Sp{n)/I, respectively 
S{n)/X which we also call basic sums. In the light of Theorem 2.3.16 this is connected 
to the number of Lyndon words. The number of Lyndon words of length n over an 
alphabet of length q is given by the first Witt formula [43, 44, 34]: 



(2.15) 



d\n 



where 



li{n) 



1 





-1)^ 



if n = 1 

if n is divided by the square of a prime 
if n is the product of s different primes. 

is the Mobius function. 

As we would like to count the number of basic sums for all sums of a given index set 
individually, this relation cannot be used: we have to count the number of Lyndon 
words belonging only to this index set. The respective number has been given in the 
same paper as the second Witt formula. 



Zn(ni, 



1 



n 



(^)!---(- 



^)!' 



n 



k=l 



(2.16) 



here denotes the multiplicity of the indices that appear in the index set. For more 
details on these aspects we refer to [34]. In other words, by Theorem 2.3.16 and with 

Sp{n)/Ip ^ i? (N) ^ Spin)/ ~ 

the number of algebraic independent sums in Sp{n)/Xp or in Spin)/ ~ (considered as 
sequences) for a given index set is 



1 (-)! 



n 



k=l 



rik 



(2.17) 
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Number of 


Weight 


Sums 


a-basic sums 


Sums -1 { 


-1} 


a-basic sums 


1 


2 


2 




1 


1 


2 


6 


3 




3 


2 


3 


18 


8 




7 


4 


4 


54 


18 




17 


7 


5 


162 


48 




41 


16 


6 


486 


116 




99 


30 


7 


1458 


312 




239 


68 


8 


4374 


810 




577 


140 



Table 2.1: Number of harmonic sums and number of sums which do not contain the 
index {—1} in dependence on their weight; respectively, the numbers of basic 
sums (a-basic sums) by which all sums can be expressed using the algebraic 
relations; compare also [5, 8]. 

again denotes the multiplicity of the indices that appear in the index set. Moreover, 
the number of algebraic independent sums in 

S{n)/I^ R{Z/{0}) 

for a given index set is 

1 (-)' * 

lnini,...,ng) = -I^A^W .niv '(^)\ ' n = J2'^k- (2.18) 

d\n ^ d. ^' ^ d, )■ k=l 

As pointed out, e.g., in [5, 6] the number of algebraic independent sums can be sum- 
marized in Table 2.1. 

Remark 2.5.1. In the difference field setting of U.T,-fields one can verify algebraic 
independence of sums algorithmically for a particular given finite set of sums; see [38]. 
We could easily verify up to weight 7 that the figures in Table 2. 1 are correct, interpreting 
the objects in S{n)/ ~ ,ie., as sequences. Nevertheless, unless we do not have a 
rigorous proof for S{n) /X = S{n)/ ^, we can only assume that the figures in Table 2.1 
give an upper bound. 



2.6 Deriving of Relations between Harmonic Sums 

In order to express the harmonic sums of a given depth in terms of a minimal set of 
harmonic sums, we will determine all algebraic relations of the harmonic sums up to 
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weight 7 following the ideas from [5]. We define the shuffle product on harmonic sums 
in the obvious way: 

Savin) luSbwin) = Sa(vUJbw){^^) + 'S'^ 

aiiLU w) 

(n), (2-19) 

for all a,b E A and v,wE:A*. 
As an example we state 

^ai,a2{n) ^ Sa^^a^in) = "5*01,02,03,04 (^) ~l" 'S'01,03,02,04 (^) ~l" 'S'oi ,03,04,02 C'^) ~l~ 

'S'o3,a4,ai,a2 l'^) + 'S'a3, 01,04, 02 l*^) + 'S'a3,ai,02,04('^)- 

We get the following relations listed up to depth 4, for further relations up to depth 6 
see [5]: 

Depth 2 

= Sai (n) LU S02 (n) - S'aiAas W " Sai (n)S'a2 (n) 



Depth 3 



= Sai{n)LuSa2,a3{n) - Sai Aa2,a3 {^) Sa2,aiAa3 (n) - Sai(n)5a2,a3(n) 



Depth 4 



= 



= 



Sai{n)inSa2,a3,a4in) - (n)- »5o2,ai Aa3,a4 in)- »S'a2,03,aiAa4 (n) 

{n)Sa2 (n) 

(n)m (n) + SaiAas ,a2 Aa4 

A'X3 ,02 )'^4 

(n) 

'S'ai Aa3,a4,a2 

(n)- »5'ai,02Aa3,a4 in) - Sa, ,03,a2 Aa4 (n)- '^03 ,ai Ao4,a2 

(n) 

^^03,01,02 Aa4 (n) 



As observed in [11] and worked out in detail in [5] these relations in combination with 
the Witt formula enable us to hunt systematically for all algebraic relations in iS(n) /X. 



Example 2.6.1. We now consider the case of harmonic sums of depth 3 as an example. 
We start with 3 different indices. There exist 6 different sums of this type: 

Sai,a2,a3{n), Sai,a3,a2{n)i 'S'o2,ai,a3 ('^)) 'S'a2,a3,ai (^)) 'S'ag, 01,02 {n), Sa3,a2,ai{n). 
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Now we use the relation we obtained for sums of depth 3 and apply it for all 6 permu- 
tations of the 3 indices. This leads to the following system of equations: 





0.2 ^^3 (^'') ^ 


ai 


a.3 ('^) + Sa2 


03 




= Sal(n)S'a2,a3('''i) + S'a2,al Aa3 (^"0 ^ 


- SalAa2,a3{'n) 


Sai 


<23 


02 (^) + 


ai 


02 ('^) "1" *S'a3 


<22 


ai(n) 


= S'al(n)Sa3,a2(?'i) + 5'a3,al Aa2 (") - 


- 5'alAa3,a2(?l) 


Sai 


<22 


a3 {^) ~\~ Sai 


a-3 


<22 ('^) "1" *S'a2 


ai 


a3(»l) 


= S'a2(n)S'al,a3(n) + Sal ,a2Aa:i{n) - 


- SalAa2,a:i{n) 


Sa2 


03 


a\ (^) "1~ Sa^ 


ai 


02 ('^) "I" Sa^ 


02 


ai(n) 


= Sa2{n)Sa3,al{n) + 5'a3,<ilAa2 (n) + 5'a2A<i3,al («) 


Sai 


12 


as iri) + Sa\ 


0,3 


12 (^) + Sas 


ai 


asW 


= 503(71)501, a2(n) + Sal,a2Aa3{n) + S'alAa3,o2 (n) 


Sa2 


ai 


a-S C*^) "1" "^12 




a\ {ri) -\- Sas 


02 




= Sa3{n)Sa2,al{n) + /So2,alAo3(n) + 'S'a2Aa3,ol (fl) 



which is equivalent to: 

/ 1 1 1 \ 

10 11 

1110 

1 1 1 

110 10 

Vooiioi/ 



^ »5ai ,a2 ,a3 (^'') ^ 
•S^ai ,a3,a2 (^) 



( Sai{n)S 

Sai i.^)Sa^,a2 i'^) H" Sa^,aiAa2 (^) ^'aj Aa3,a2 (j^) 
Sa2 (^)'S^ai ,a3 (^) "I" ,a2Aa3 ('^) H" SaiAa2,a3 C?^) 

*5'a.2 (''^^) *5'a3 ,ai ('^^) "h 5'a3,aiAa2 

(n) + S'a2Aa3.ai 

ai ,a2 

(n) + 5ai,a2Aa3 + Aa3 ,02 ("■) 

V Sophs' a2,ai(7i) + 5*02,01 Aa3 (ra) + 5*02 Aa3 ,01 (jl) / 



Since the rank of this linear system is 4, all sums with 3 different indices can he expressed 
by two chosen sums of depth 3 and sums of lower depth. To be more precise, solving 
this system leads to 



^oi A03,02 {n) + So3 ,ai Aa2 

(n) + Sai(n)So3,a2 W - »S'a3,a2,ai in), 

(n) Sai Aa3,a2 

(n) - S'o3(n)Soi,a2(«) + 5o2(fl)>S'ai,03(n) +S'o3 ,ai ,a2 in), 



~\~Sa3,aiAa2 

(n) + 5oi(n)S'o3,a2 (n) Sa3,ai,a2 

(n) Sas,a2,ai (n) 

For two different indices we get the following system of equations: 



Sai 


12 


03 (n) 


— Sai,a2Aa3 


Sai 


as 


a2H 


— SaiAa2,o.2 


Sa2 


oi 


03 H 


— SaiAa2,o.3 


Sa2 


"3 




— SaiAa2,a2 




^ Sai ,0.1 ,a2 (^) 
Sai,a2,ai (^^) 
\ *^a2,oi ,ai (^) 



/ Sai{n)Sai,a2{n) + Sai ai Aa2 

(n) + SaiAai,a2{n) 
Sai {n)Sa2,ai {n) + Sa2 ai Aoi {n) + 5'aiA02,ai(n) 
\ Sa2{n)Sai,ai{n) + Sai ,ai Aa2 

(n) + 5oiAo2,ai(n). 



Since the rank of this linear system is 2, all sums with 2 different indices of depth 3 
can be expressed by one chosen sum of depth 3 and sums of lower depth: 



Sa 



" [Sai Aai,a2 (^) SaiAa2,o.i (j^^ ~^ Sai,ai Aa2 (^) H" Sai {_^^Sai ,0,2 (^) 



Sa2 ,ai Aai 

(n) - Sai {n)Sa2,ai («) + S'02 ,ai ,ai (n)] + Sa2 ,ai ,ai in), 

Sai,a2,ai in) = 

iJaiAa2 ,ai oiAoi {n) + Sai{n)Sa2,aiin)-2Sa2 ,ai ,ai in). 

For a sum of depth three with 3 equal indices we obtain: 

1, 



(2.20) 



Sai ,o,i ,ai (^) 



, [*^aiAai,ai (^) H~ 'S'a,-|^ (n) + S'oi(n)S'oi,oi(n)l- 



Remark 2.6.2. A different method which can be used for equal indices and is described 
in [5] leads to: 



1 3 

~ ^ \Sai (^) H~ ^Sai * SaiAai 

(n) + 2 * Sa 1 Aai Aai (n)]. 
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Applying the second Witt formula (2.16), i.e., 

Zn(l,l,l) =2 
^n(2,l) = 1 
^n(3) = 

we can see that we found the maximal number of algebraic independent sums. In other 
words, this strategy solved the problem of weight 3 completely. 



Using HarmonicSums , these relations can be computed as follows: 



HarmonicSums session. 

in[5]:= GetDependentSSums[{l, 1, 1}] 

sums: 6 

dependent sums: 4 

Out[5l= {S[al_,a2_, a3_,n_] -> S[a3, n]S[al, a2, n] + S[al, SP[a2, a3], n] - S[al, 7i]S[a3, a2, ri] - 
S[a3, SP[al,a2],n] + S[a3, a2, al, n]/; MyMcmbcrQ[al, a2, a3, al, a2, a3]&&al / a2 / a3, 
S[al_,a3_,a2_,n_] ^ S[al, n]S[a3, a2, n] + S[a3, SP[al, a2], n] + S[SP[al, a3], a2, n] - 
S[a3, al, a2, n] — S[a3, a2, al, n]/; MyMemberQ[al, a2, o3, al, o2, a3]&&al ^ a2 / a3, 
S[a2_,al_,a3_,n_] -S[a3,n]S[al,a2,n] + S[a2,n]S[ol,a3,n] + S[SP[ol, a2], a3, n] - 

S[SP[al, a3], o2, n] + S[o3, ol, a2, n]/; MyMemberQ[al, o2, o3, al, a2, o3]&&al ^a2^ a3, 
S[a2_,a3_,al_,n_] -S[o2, n]S[al, o3, n] - S[al, SP[a2, o3], n] + S[al, n]S[a2, a3, n] + 

S[a2,SP[al,o3],n] + S[ol,n]S[o3,o2,n] + S[o3, SP[al, o2], n] + S[SP[ol, a3], o2, n] - 
S[a3,al,a2,n] - S[a3, o2, ol, n]/; MyMemberQ[al, a2, o3, ol, a2, a3]&&al ^ a2 ^ o3} 

in[6]— GetDependentSSums[{2, 1}] 

sums: 3 

dependent sums: 2 

Out[6]= {S[ol_,ol_,o2_,n_] ^S[al,n]S[al,a2,n] + ^S[al, SP[ol, o2], n] - ^S[al, n]S[a2, ol, n] - 
^S[a2,SP[ol,ol],n] + ^S[SP[al, ol], a2, n] - ^S[SP[al,a2],al,n] + S[o2,ol,al,n]/;ol ^ a2, 
S[al_,a2_,al_,n_] S[ol, n]S[o2, ol, n] + S[a2, SP[ol, ol], n] + S[SP[ol,o2],ol,n] - 

2S[o2,ol,ol,n]/;ol o2} 
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2.7 Application of the Relations and Algebraic 
Simplification 

There are now basically three ways to use the strategy presented in the last section. 

2.7.1 First Strategy: Fixed Tables 

We can fix basic sums and compute tables up to a certain weight. Then we can use 
these tables to reduce expressions in harmonic sums to expressions where only the 
chosen basic sums appear. Tables up to weight 6 are produced in [5] whose basic sums 
are of particular interest in particle physics. For the package Harmonic Sums these 
relations were recomputed and can now be used in the following way. 

HarmonicSums session. 

Relations up to weight 6 are included in the package. 

in[7].= expr = S[l, 2, -2, n] + S[-l, 2, 1, n] + S[-l, 2, -2, n] + S[-2, 2, 1, n] + S[-2, 1, 2, n] 
+ S[l, -2, -1, n] + S[-l, 2, 1, n] + S[2, -1, 2, n] + S[-l, 2, 2, n]; 

in[8]:= ReduceToBcisis[expr] 

n]S[— 1, n] + S[— 4, n]S[l, n] +S[— 2, ri]S[l, 'm]S[2, n]+S[—2, n]S[3, n] + 
^S[4, n] + S[-l, n] [S[2, nf + S[4, n]] - S[-4, -1, n] - S[-4, 1, n] + S[-3, -2, n] + S[-3, -1, n] + 
2S[-3,2,n] + S[l,n]S[-2,-l,n] + S[2, nlS[-2, 1, n] - S[-l,n][S[-4,n] + S[-2, n]S[2, n] - 
S[2,-2,n]] - S[l,n][S[-4,n] + S[-2,n]S[2,n] - S[2,-2,n]] + S[-2, n][S[-l, n]S[2, n] - 
S[2, -1, n]] - S[2, n]S[2, -1, n] - 2S[-2, n]S[2, 1, n] - S[4, -1, n] - S[-2, -1, 1, n] - S[-2, 1, -1, n] + 
S[-2,2,-l,n] + 2S[-2,2,l,n] - S[-l,2,2,n] + S[2,l,-2,n] + 2[S[-3,l,n] + S[2,-2,n] + 
S[-l, n]S[2, 1, n] - S[2, -1, 1, n] - S[2, 1, -1, n]] 



The disadvantage is however, that for certain applications it might be more reasonable 
to chose other basic sums. In this case we cannot use this strategy. 

2.7.2 Second Strategy: Partly Fixed Tables 

We can compute tables for all possible index sets up to a certain depth, which is done 
in [5] up to depth 6. These tables have been extended up to depth 7 during the work on 
this thesis. The computation of these tables took about one and a half week CPU time 
on a 2.3 GHz machine, and ended in a file with about 300 Mbytes. We can use these 
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tables to reduce expressions in harmonic sums to expressions where for each depth and 
index set at most that many sums appear as is given by the second Witt formula. 



HarmonicSums session. 






First we have to load the table where the relations are stored: 






in[9]- DeflneRelationTable["RelTab.m"] 






in[ioi:= expr = S[l, 2, -2, n] + S[-l, 2, 1, n] + S[-l, 2, -2, n] + S[ 


-2, 2, 


1, n] + S[-2, 1, 2, n] 


+ S[l, -2, -1, n] + S[-l, 2, 1, n] + S[2, -1, 2, n] + S[-l, 2, 


2, n] 




in[ii]:= ReduceToBasis[expr, Dynamic —i- Partly] 






Out[ii]= S[-5,n] + S[-3,n]S[2,n] + ^S[-l,n][S[2,n]^ + S[4,n]] - S[ 


-4,1, 


n] + S[l,n]S[2,-2,n] - 


S[-2, n]S[2, 1, n] + S[3, -2, n] + S[-2, 1, 2, n] + 2S[-2, 2, 1, n] + S[ 


-1,2, 


-2,n] + 2S[-l,2,l,n] + 


S[l,-2,-l,n]-S[2,2,-l,n] 







However there is the disadvantage that although for each depth and index set the 
relations can be calculated quite generally, we have to fix in advance certain sums 
which have to be in the basis. We can easily see this in the following example: 
We again consider sums of depth 3 with 2 different indices, so we look at the sums 
'S'ai, 01,02 ,02,01 1'^)) "^02,01,01 ('^)- According to the W^itt formula there should be 

just one of the sums in the basis. In (2.20) we took Sa2,ai,aiin) to express the other 
sums. Therefore for example the sum S'i,i,2(w) will always be replaced by 52,1,1 (n) and 
sums of lower depth. Sometimes however it could be more convenient to take S'2,i,i(n) 
as a basis element. 

2.7.3 Third Strategy: Dynamic Reduction 

If we look at an expression consisting out of several different harmonic sums, we could 
first look at all the different index sets of harmonic sums that appear in the expression. 
For each index set J we can use the Witt formula to get the number n j of basis elements, 
i.e., the basic sums. Then for each index set we can use our method to set up a system 
of equations where we choose exactly n of the permutations (and therefore n harmonic 
sums) to be in the basis. Here we will try, as far as possible to choose sums which are 
already in the expression, since we do not want to introduce new sums. Sometimes 
however we will have to choose also sums which are not in the expression. By solving 
the equation system the rest of the sums can be expressed by harmonic sums out of the 
basis and sums of lower depth. So we end up with an expression, where for each index 
set J at most nj different multiple harmonic sums appear. Again we can illustrate this 
by an example. 
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Example 2.7.1. Consider an expression where the sums Sai^ai^a2(^^^ and Sai.a2^aii'^^ 
appear. We will use one of the sums to express the other. Since the Witt formula gives 
1, this is possible. In the reduced expression there will he either the first, the second 
sum or no sum with this index set. Using, e.g., the second strategy we would use (as 
defined in Example 2.6.1) Sa^.a^aiin), that is not in the original expression. 



HarmonicSums session. 






Of course, now we do not have to load any tables, since they are computed online. 






in[i2]:= expr = S[l, 2, -2, n] + S[-l, 2, 1, n] + S[-l, 2, -2, n] + S[-2, 2, 1, n] + S[ 


-2, 1, 


2, n] 


+ S[l, -2, -1, n] + S[-l, 2, 1, n] + S[2, -1, 2, n] + S[-l, 2, 2, n]; 






in[i3]:= ReduceToBasis[expr, n, Dynamic — > True] 






Out[i3]= S[-5, n]+S[-3, n]S[2, n] + ^S[-l, n] [S[2, nf+S[4, n]]+S[l, n]S[-2, 2, n]+S[-2, 3, n]- 


-S[l,- 


-4,n]- 


S[-2, n]S[l, 2, n] - ^S[2, -3, n] - ^S[2, n]S[2, -1, n] - ^S[4, -1, n]+S[-2, 1, 2, n] +S[- 


1,2,- 


-2,n] + 


2S[-1, 2, 1, n] + S[l, -2, -1, n] + 2S[1, 2, -2, n] + ^S[2, -1, 2, n] 







This strategy works quite fine (which means fast) for index sets where the number of 
permutations is not too high. For example there are 60 permutations of an index set 
of the type {01,01,02,03,04} in such a situation the strategy works really fast. But 
if we consider the case of an index set of the type {oi, oi, 02, 03, 04, 05, ag} , there are 
2520 permutations. If we have to set up and solve the system of equations in this case, 
this will take hours. 



2.7.4 Algebraic Simplification 

We use the notions of [13] concerning canonical algebraic simplification. Especially 
the first strategy can be used to obtain equivalent but simpler objects and to compute 
unique representations for equivalent objects in S{n)/I; here I is again the ideal in 
(2.14). For s,t E S{n), we say s < t if the number of multiple harmonic sums in s 
which arc not basic sums is smaller or equal than the number of non basic sums in t. 
Now consider the following problem: 

Given: s G 5(n). 

Find: t G S{n) such that s + I = t + I and such that for all a G S{n) with a + I = t + I: 
t < a. 
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Now let T be the procedure which uses the first strategy. T maps S{n) to S{n). For 
ah s,t € S{n) with s + I = t + I wc have the equahty T{s) = T{t) in S{n) by 
construction; T{s) / T(i) would mean that algebraic relations are missing. Hence we 
can check equivalence easily. 

Additionally, the number of non basic harmonic sums in T{s) is zero; hence the above 
problem of simplification is solved completely. 



2.8 Relations between Similar Structures 

In the same way one can find relations for Euler-Zagier sums, here the product (2.4) is 
the quasi shuffle product discussed in [21]: 

e* w = w * € = w, for all w E A* 
au*bv = a{u* v) + b{u* v) + [a,b]{u* v), fox all a,b e A;u,v & A* 

with [a,b] = a Ab for all a, 6 G Z. 

Again the number of Lyndon words gives the number of algebraically independent 
Euler-Zagier sums viewed as elements in a ring similarly as S{n)/I for harmonic sums 
and it gives at least an upper bound for the number of algebraic independent Euler- 
Zagier sums viewed as sequences. In the case of just positive indices it was proven 
in [27] that the number of Lyndon words is not only an upper bound for the number 
of algebraically independent Euler-Zagier sums viewed as sequences, it is in fact the 
number of algebraically independent Euler-Zagier sums viewed as sequences. It is likely 
that the same holds for index sets where we also allow negative indices. Like for har- 
monic sums the package HarmonicSums provides a procedure to calculate the relations 
automatically: 

HarmonicSums session. 

in[i4]:= GetDependentZSums[l, 1, 1] 

sums: 6 

dependent sums: 4 

Out[i4]= {Z[al_,a2_,a3_,n_] Z[a3,n]Z[al,a2,n] - Z[al,SP[a2,a3],n] - Z[al,n]Z[a3,a2,n] + 

Z[a3,SP[al,a2],n] -FZ[a3,a2,al,n]/;MyMemberQ[al,a2,a3,al,o2,a3]&&al ^ a2 ^ a3, 
Z[al_,a3_,o2_,n_] Z[ol,n]Z[a3,o2,n] - Z[a3, SP[al, o2], n] - Z[SP[al,o3],a2,n] - 

Z[a3, ol, a2, n] - Z[o3, a2, al, n]/; MyMemberQ[ol, o2, a3, al, a2, o3]&&:ol a2 o3, 
Z[a2_,al_,o3_,n_] -Z[a3,n]Z[al,o2,n] -|- Z[a2,n]Z[al,o3,n] - Z[SP[al,o2],a3,n] -|- 

Z[SP[ol, o3], a2, n] + Z[a3, al, o2, n]/; MyMemberQ[al, a2, o3, ol, o2, a3]&&al ^ a2 ^ a3, 
Z[a2_,a3_,ol_,n_] -Z[a2,n]Z[al,o3,n] -|- Z[al,SP[o2,a3],n] + Z[ol,n]Z[o2,a3,n] - 

Z[a2,SP[ol,o3],n] + Z[al, n]Z[o3, o2, n] - Z[o3, SP[al, a2], n] - Z[SP[al, a3], o2, n] - 



Algebraic Relations between Multiple Harmonic Sums 



25 



Z[a3,ol,o2,n] - Z[o3,o2,al,n]/;MyMemberQ[ol,o2,o3,al,a2,o3]&&al ^ a2 ^ o3} 



We can generalize the harmonic sums to S-sums [28]. For Oj G Z and G M we define 

Sai,...,ak;xi,...,Xk{l^) = ^ (2-21) 

The product of two S-sums yields: 



— X X 

n>ii>i2>--->ife>l ^ ^ 



Sai,...,ak;xi,...,Xk (n) * S'6i,...,bj;yi,...,2^j (n) — ^ 'S'a2,->afc;x2,...,a;fc (0 * Si,^,...,bi;yi,-,yi (*) 

n 



1=1 

~ ^ ,-oi+6i '^a2,.--,afc;a;2,.--,^fc(*) * ^b2,-,bi;y2,-,yi{ 

i=l ^ 

(2. 

Again recursive application of (2.22) leads to a linear combination of single S-sums. 
As showed in [28], ideas from the previous sections can be carried over as follows. We 
consider an alphabet A, where pairs (m, x) form the letters. Moreover, we define the 
addition and the multiplication of words by 

[(mi, xi), (m2, X2)] = (mi -|- m2, xi * X2) 

and by 

(mi, xi)it;i * (m2, ^2)^2 = {mi,xi){wi * (m2, X2)w2) + (m2, X2){wi * (m2, X2)w2) 

-[(mi, xi), (m2, X2)\{wi * W2). 

Note that this multiplication is identical to the multiplication of S-sums (2.22). The 
degree of a letter is |(m,a;)| = m. Since (S0)-(S3) hold for [•,•], the S-sums form a 
quasi-shuffle algebra. Finally, we emphasize that we implemented algorithms to find 
algebraic relations between S-sums following the ideas presented in Section 2.6. 



HarmonicSums session. 

in[i5]— GetDependentGeneralizedSSums[2, 1] 

sums: 3 

dependent sums: 2 

Out[i5]= {S[al_, al_, a2_, a;l_, a;l_, a;2_, n_] — >■ ^S[al,xl,n]Q[al,a2,xl,x2,n\ + ^S[al,al + 
a2,a;l,a;la;2,n] + ^S[2al,a2,a;l^,a;2,n] — ^S[ol,a;l,n]S[o2, al,a;2,a;l,n]— ^S[a2, 2al,a;2,a;l^,n] — 
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^S[al + a2,al,xlx2,xl,n] + S[a2, al, ol, a;2, xl, xl, n], 

S[al_, a2_, al_, a;l_, x2_, a;l_, n_] — > S[al, a;l, n]S[a2, al, a;2, xl,n] + S[o2, 2al, a;2, xl'^, n] + S[al + 
a2, al, xlx2, xl, n\ — 2S[a2, ol, ol, x2, xl, xl, n]} 



In a similar way we can generalize Euler-Zagier sums to Z-sums [28]: 

Zai,...,ak;xi,...,xu{n) = ^ ' ' ' ^ ■ (2.23) 

n>ii>i2>->2fc>0 ^ 

The product of two Z-sums yields: 

n j 

Zai,...,ak;xi,...,Xk{n) * ■^6i,...,bi;i/i,...,2/i (?^) = ^ ^■^a2,.--,afe;a:2,.--,a;fc (*) * -^6i,...,6;;2/i,...,2/i " 1) 



(ai * 61)' 



«=1 



1=1 



Similar as for harmonic sums 
and Z-sums [28]. 


and Euler-Zag 


ier sums we can convert between S-sums 


HarmonicSums session. 






In[16]:= ZToS[Z[l, 3, 4, {2, 1, 1}, 


n]] 




Out[i6]= S[8, {2}, n] - S[l, 7, {2, 1}, n 


-S[4,4,{2,l},n 


+ S[l,3,4,{2,l,l},n] 


ln[17]:= SToZ[S[l, 3, 4, {2, 1, 1}, 


n]] 




Out[i7i= Z[8,{2},n] + Z[l,7,{2,l},n] 


+ Z[4,4,{2, l},n] 


+ Z[l,3,4,{2,l,l},n] 
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Chapter 3 

Harmonic Polylogarithms and the 
Differentiation of Harmonic Sums 

The bigger part of this chapter will deal with harmonic polylogarithms. Harmonic poly- 
logarithms were first introduced in [33] and are covered by Poincare iterated integrals 
[31, 24]. For the definition and basic properties of the harmonic polylogarithms we will 
follow more or less [33] , however sometimes we will go more into detail (for example we 
will provide detailed proofs), and we will skip aspects which are not important for our 
considerations. 

We discuss harmonic polylogarithms since it will turn out that they are connected to 
multiple harmonic sums via a special extension of the Mellin transform, which is fre- 
quently used in particle physics. The Mellin transforms of harmonic polylogarithms are 
expressions in multiple harmonic sums, while the inverse Mellin transforms of multi- 
ple harmonic sums are expressions in harmonic polylogarithms. Hence we can use the 
Mellin transform to convert between multiple harmonic sums and harmonic polyloga- 
rithms. As a direct consequence we can use the inverse Mellin transform to construct 
analytic continuations of multiple harmonic sums. An algorithm to compute the Mellin 
transform of harmonic polylogarithms is given in [33]; we will analyze this algorithm 
in detail and will verify its correctness. Similarly, an algorithm for the computation 
of the inverse Mellin transform of multiple harmonic sums is given in [33]; again we 
will work out the correctness of this algorithm and give additional insight from a more 
computer algebra point of view. In order to execute these algorithms, we will also need 
to consider harmonic sums at infinity and harmonic polylogarithms at one. Both do 
not have to be finite, but it will turn out that they are connected via power series 
expansions. 

These algorithms are originally implemented in Vermaseren's package harmpol in the 
computer algebra system form [33]. Harmpol also provides procedures to manipulate 
with harmonic polylogarithms. In addition, there is the Mathematica-package hpl [26], 
in which basic procedures for manipulating harmonic polylogarithms are available. 
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inverse Mellin transform 



ellin transfo: 







i 

8 




i 



power series expansion 



Figure 3.1: Connection between harmonic sums and harmonic polylogarithms. 

Since we can use the inverse Mellin transform to compute analytic continuations of 
multiple harmonic sums, we can define a differentiation on multiple harmonic sums as 
done in [6, 8, 7]. 

In the last part of this chapter we give an algorithm to do this differentiation auto- 
matically. In the package HarmonicSums the Mellin-transform and the inverse Mellin- 
transform together with procedures for harmonic polylogarithms which are also needed 
are implemented. In addition, the differentiation of multiple harmonic sums is imple- 
mented. Along with this differentiation new relations between multiple harmonic sums 
are going to emerge, and these relations will reduce the algebraic basis computed in the 
previous chapter. 



3.1 Definition of Harmonic Polylogarithms 



In order to define harmonic polylogarithms as introduced in [33], we first define the 
function / : {0, 1, -1} x (0, 1) M by 



f{0,x) 



1 



fihx) 



X 

1 



1 



\ + x' 
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Harmonic polylogarithms arc fmictions in C°^((0, 1)) and are defined inductively: Let 
m = {rriw, m^^i, . . . , rrii), rrii G {0, 1, —1} be a vector of length t/; G N and let 0^ be 
the vector of length w whose w components are all equal 0. For x G (0, 1) we define 
the harmonic polylogarithm H(m, x) : 

H((),x) = 1, 
H(0„,x) = ^(loga^r, 

B.{m,x) = / f {mw,y)'ii{{m.yj-i,mw-2,---, mi), y)dy if rriy^ Oyj- 
Jo 

The length of the vector m is called the weight of the polylogarithm H(m, x). 
Examples are 

H 



{{i),x) = r^=-iog(i-x), 

Jo ^ - y 

H((-l),x) = r^=log(l + x), 
Jo J- "T 2/ 



l + y Jo 1 + 2/ 



dy. 



Remark 3.1.1. (compare [33]) It follows from the definition that ifm ^ Ow, H(m, 0) = 
0. If my, ^ 1 or if mw = 1 and = for all v with v < w then H(m, 1) is finite. In 
the remaining, cases, i.e., my, = 1 and my ^ for some v with v < w, Iim^_^i- H(m, x) 
behaves as a combination of powers o/log(l — x) as we will see later in more detail. 
We define H(m, 0) := Iim^_^o+ H(m, x) and H(m, 1) := Iimj._>.i- H(to, x) if the limits 
exist. 

Notation 3.1.2. Prom now on we will write iimvj,m.u,-i,-,mi{x) for 
H((m^,m^_i, . . .,mi),x). 

Remark 3.1.3. For the derivatives we have for all x G (0, 1) that 

— Hm(a^) = f{my,,x)}lm^-i,my:-2,-.-,mi{x). 



Subsequently, we present some identities between harmonic polylogarithms of the same 
argument. These identities will lead to a formula, which we can use to rewrite the 
product of harmonic polylogarithms to a sum of single harmonic polylogarithms. 

Lemma 3.1.4. For x G (0, 1), q eN and mi G {0, 1, -1} , 

~ ^^im2,mi{x)iim3,...,mq{x) 
+ ^m3,m2,mi{x)^m4,...,mq{x) 

- ...-{-i)mmq,...,mAx). (3.1) 
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Proof. Case 1, (mi, . . . , rriy) = Oq. 
From the definition we get 

^mu...,m^{x) = ^logx^ 
and the right hand side of (3.1) yields 

^°g^ (g^l)! x-^^ log2 x ^^^^^, log«-2 X 
log^ ^^A^ log""^ a; (-1)^^ log'^ X 

1=1 

Case 2, (mi, . . . , m^) 7^ Oq. 

Using the definition and partial integration we get 

PX 

Hmi,...,m,(a;) = / /(mi,y)H 

1712, ■■■,mq {y)dy 

Jo 

= H^i(x)H^2,...,m,(a;) - / f{m2,y)B.^3,...,mq{y)dy. 

Jo 

Continuing by partial integration we finally get the right hand side. □ 



We can use this lemma to get the following identity. 

Lemma 3.1.5. For x G (0, 1), a G {0, 1, -1}, q eN and rrii G {0, 1, -1}, 

Ha(2^)Ilmq,...,mi (•^) — Ha,mq,mq_i,...,mi (•^) 

"I" Hmq,a,mq_i,...,mi (^) 

+ ^mq,mq-\,a,mq-2,...,iniix) 

+ ••• 

+ ^mq,mq-i,...,mi,a{x). (3-2) 



Proof. Let q = \: By (3.1) we get Ha(a;)H^j (x) = 'B.a,mAx)- Assume (3.2) holds for q 

^a{x)^mq,...,m\{x) = ila,mq,mq-i,...,mi{x) 

+ ^mq,a,mq-i,...,mi{x) 

+ ^mq,mq-i,a,mq-2,—,mi{x) 

+ ■■■ 

+ ^mq,mq-i,...,mi,a{x). 



Harmonic Polylogarithms and the Differentiation of Harmonic Sums 



31 



Multiply both sides by f{mq^i,y) and integrate: 

/("ig+i, y)Ha(2/)H^„...,mi {y)dy 

{y)dy 

Jo 

{y)dy 

Jo 

{y)dy 

Jo 

{y)dy- 

Jo 

By partial integration on the left hand side and the definition of harmonic polyloga- 
rithms we get: 

/•X 

Ra{x)Rrnq,...,mi{x) - / fia,y)Rmq+u-mi{y)dy 

Jo 

— R-mq+i,a,mq,...,mi{x) + Hmq+i,mq,a,mq-i,...,mi{x) 
''rR-mq+i,mq,mq-i,a,mq-2,-..,mi{x) + • • • + iimq+i,...,mi,a{x)- 

Hence (3.2) holds for q + 1. □ 



Similar to multiple harmonic sums, the shuffle product on words is extended to harmonic 
polylogarithms. 

Definition 3.1.6. Let a,b£A*. If 

am 6 = cidi H \-Cndn 

for di E A* , Ci E M. then we define 

HaLU6(a;) := ciHrf,(x) H l-CnHd^(x). 

Definition 3.1.7. Let x G (0, 1), a = (ai, 02, . . . , Up) and b = (61, 62, . . . , bq), where 
ai, hi G {0, 1, —1}, p,q € N. Let us interpret a and b as words, so a = aia2 • • • ap and 
b = bib2 ■ ■ ■ bq. We define the shuffle product of two harmonic polylogarithms by 

H„(x)LuHfe(x) := {a2---apLU b){x) + ^bi{ainb2-bq){x) 

where }lcic2-cn{x) ■= R-ci,c2,:;Cnix) for ci G {0,1,-1} and where lu is defined as in 
Definition 2.3.4. 
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Examnle 3 18 For r F (0 1) 




Hi,o-i(a::)LuHo,i(x) = 


Ho,i,o -i,i(^) + Ho,i,o,i -i(^) 


+ 


2Ho,i,i,o,-i(.'r) + Hi,o -i,o,i(a;) 


+ 


2Hi,o,o + 2Hi,o,o,i 


+ 


Hi,o,i,o,-i(a;)- 


HarmonicSums provides a procedure to carry out this product: 


HarmonicSums session. 




in[i8]:= ProductH[H[l,0,-l,x],H[0,l,x]] 




Out[i8l= H[0, 1, 0, -1, 1, x] + H[0, 1, 0, 1, -1, x] 


+ 2 H[0, 1, 1, 0, -1, x] + H[l, 0, -1, 0, 1, x] 


+ 2 H[l. 0. 0, -L 1, x] + 2 H[l, 0. 0. 1, 


-1, x] + H[l, 0, 1. 0. -1. x] 



Theorem 3.1.9. Let x G (0, 1), a = (op, Up-i, . . . , oi) and b = {bq, bg^i, . . . , bi), where 
ai, bi G {0, 1, —1}, p,q EN. For the product of two harmonic polylogarithms we get 

H„(x)Hb(x) = H„(x) mH(.(a;). (3.3) 
Proof. We aheady showed the case for p = 1. Now assume (3.3) holds for p > 1 : 

^ap,ap-i,...,ai (a^)H^,6,_i,...,bi (x) = iiap,ap-i,...,ai (a;) LU H^_fe^_-^^...^bj (x). 

Multiply both sides by /(ap+i,x) and integrate: 

X 

f{ap+l,y)^ap,ap-i,...,ai{y)^b^A-U...Miy)dy 

fiap+i,y) (H 

ap,ap-i,...,ai 

Partial integration of the left hand side gives 

Hap+i,ap,...,ai {x)ilbq,bq-i,...,bi {x) 
px 

(y)^H6^_i,b,_2,...,6i(2/))f^y- 

JO 

Let us abbreviate a = Cp+iOp • • • ai, hence we may write 

Ha(x)mH6(a;) = Ha^+^,ap,...,aAx)^^b„b,-^,...,bAx)- 

Expanding the first shuffle product, using the definition of the harmonic polylogarithms 
and moving it to the right hand side gives a sum of all those harmonic polylogarithms 
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which appear in Ha(a:) luH5(,t) and start with bq, while expanding the second shuffle 
product gives a sum of all those harmonic polylogarithms which appear in Ha(x) LuH5(a;) 
and start with a^+i. Since all harmonic polylogarithms in Ha(a;) LuHft(a;) either start 
with bq or Op+i we completed the proof. □ 



The following remarks are in place. 

Remark 3.1.10. For harmonic polylogarithms which are finite at x = and x = 1, 
(3.3) can be extended to x £ [0, 1]; however if one of the harmonic polylogarithms is 
not finite, (3.3) does not hold (compare [33]). 

Remark 3.1.11. The product of two harmonic polylogarithms of weights wi and W2 
can be expressed as a linear combination of {wi+W2y./{wi\w2\) polylogarithms of weight 

W = Wi + W2. 

Remark 3.1.12. We say that a harmonic polylogarithm tlai,a2,:.,ap{x) has trailing zeros 
if there is an i E {1, 2, . . . ,p} such that for all j G N with i < j < p, aj = 0. 
Likewise, we say that a harmonic polylogarithm tlai,a2,.-,apix) has leading ones if there 
is an i E {1, 2, . . . ,p} such that for all j G N with 1 < j < i, aj = 1. 



We can use (3.2) now to single out terms of logx from harmonic polylogarithms whose 
indices have trailing zeros. Let a = in (3.2); then after using the definition of harmonic 
polylogarithms we get: 

log{x)Rrng,...,mAx) = Ho,m„...,mi (a::) + H^„0,m,_i,...,mi (a;) 
+Hmq,mq_i,0,mq_2,—,"ii (•^) + ■ ■ ■ + Hmg,...,mi,o(a^)) 

or equivalently: 

Hm,,...,mi,o(a^) = log(x) * H,rig,...,mAx) " ^0,mq,...,mAx) 
~^mq,0,niq-i,...,mi{x) — ■ ■ ■ — amq,...,0,mi{x) ■ 

If mi is as well, we can move the last term to the left and can divide by two. This 
leads to 

iimq,...,Ofi{x) = \ (log(a;) * ^mq,...,m2,Q{x) - Ho,,n,,...,m2,o(a;) 
~H^^_0,mq_i,...,m2,o(a;) — • • • — H^^^..._o,m2(a^)) ■ 

Now we can use (3.2) for all the other terms, and we get an identity which extracts 
the logarithmic singularities due to two trailing zeros. We can repeat this strategy as 
often as needed in order to extract all the powers of log(a;) or equivalently Ho(a;) from 
a harmonic polylogarithm. 
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Example 3.1.13. For x G (0, 1), 

1 



Hi-i,o,o(a;) = ;,Ho(x) Hi - Ho(x)Ho,i 



2 

- Ho(x)Hi,o -i(x) + Ho,o,i -i(a;) + Ho,i,o 
+ Hi,o,o-i(a;) 

= ^ log^(x)Hi _i(x) - log(x)Ho,i _i(x) 

- log(x)Hi,o -i(x) + Ho,o,i,-i(a;) + Ho,i,o,-i(a;) 
+ Hi,o,o-i(a;). 



This can be done with HarmonicSums as follows. 



HarmonicSums session. 




in[i9]:= RemoveTrailingO[H[l, -1, 0, 


X]] 


Out[19]= — - B[0,x 


H[0,l,-l,a;] - H[0,x]H[l,0,-l,x] + H[0, 0, 1, -1, x] + 


H[0, 1, 0, -1, x] + H[l, 0, 0, -1, x] 





Remark 3.1.14. In our implementation we prefer to use H[0,x] instead of Log[x]. 

Remark 3.1.15. We can decompose a harmonic polylogarithm iimq,mq-i,...,mi{x) in a 
univariate polynomial in Ho(x) with coefficients in the harmonic polylogarithms without 
trailing zeros. If the harmonic polylogarithm has exactly r trailing zeros, the highest 
power o/Ho(x), which will appear, is r. 



Similarly, we can use (3.2) to extract powers of log(l — x), or equivalently Hi(x) from 
harmonic polylogarithms whose indices have leading ones and hence are singular around 
X = 1. Let a = 1 in (3.2); then after using the definition of harmonic polylogarithms 
we get: 

log(l - x) * Bmq,...,miix) = Hi,^^,...,^j(x) + Ilmq,l,mq-u-,miix) 
''r^mq,mq-i,l,mq-2,—,mi{x) + • • • + iimq,...,mi,l{x), 

or equivalently: 

iil,mq,...,mAx) = log(l - x) * B.mq,...,rmix) - H^<„l,m,_i,...,mi (a?) 
~^imq,mq-i,l,...,mi{x) — ■ ■ ■ — iimq,...,mi,l{x) ■ (3-4) 
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If niq is 1 as well, we can move the second term to the left, and divide by two. This 
leads to 

Hl,l,...,mi(a;) = ^ (log(l -x)* Hi,^^_i,...,^i(x) - til,m,-ul,..;mi{x) 
— ■ ■ ■ — tli^rnq-i,...,mi,l{x)) . 

Now we can use (3.2) for all the other terms and we get an identity which extracts the 
singularities due to two leading ones. We can repeat this strategy as often as needed 
in order to extract all the powers of log(x — 1) or equivalently Hi(x) from a harmonic 
polylogarithm. 

Example 3.1.16. For x G (0, 1), 

Hi,i,o,-i(ic) = -Ho-i(a;)Hi,i(a;) - Hi(x)Ho 

- Hi(x)Ho,i,_i(x) + Ho,-i,i,i(x) 

+ Ho,i -i,i(x) + Ho,i,i -i(x) 

= * Ho,-i(x)(- log{x - + log(a; - l)Ho,-i,i(x) 

+ log(a; - l)Ho,i + Hq -i,i,i(x) 



+ Ho 


1,-1,1 


(x) + Ho,i,i,-i(x). 


HarmonicSums session. 






in[2o]:= RemoveLeadingl[H[l, 1, 


0, -1 


X]] 


Out[20]= H[0, — l,x]H[l, l,a;] — 


H[l,a; 


H[0,-l,l,a;] - H[l,a;]H[0,l,-l,a;] + H[0, -1, 1, 1, a;] + 


H[0, 1,-1,1, a;] + H[0, 1,1,- 


-l,x] 





Remark 3.1.17. As before we prefer to use H[l,x] instead of Log[l-x] in our imple- 
mentation. 

Remark 3.1.18. We can decompose a harmonic polylogarithm Hmg,m,_i,...,mi(a;) in a 
univariate polynomial in Hi (x) with coefficients in the harmonic polylogarithms without 
leading ones. If the harmonic polylogarithm has exactly r leading ones, the highest 
power o/Hi(x), which will appear, is r. So all divergences of harmonic polylogarithms 
for X — )• 1 can be traced back to the basic divergence o/Hi(x) = — log(l — x) for x — >• 1. 

Remark 3.1.19. We can combine these two strategies to extract both, leading ones and 
trialling zeros. Hence, it is always possible to express a harmonic polylogarithm in a 
bivariate polynomial in Ho(a;) = log(a;) and Hi(a;) = log(l — x) with coefficients in the 
harmonic polylogarithms without leading ones or trailing zeros, which are continuous 
on [0, 1] and finite at x = and x = 1. 
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Example 3.1.20. For x G (0, 1), 

Hi,_i,o(x) = -Ho(x)H_i,i(x)+Hi(x)(H_i(x)Ho(x) 

- Ho,-i(x)) +Ho,-i,i(.x) 

= — log(a;)H_i^i(a;) + log(l — a;)(H_i(x) log(a;) 

- Ho,-i(x)) + Ho,-i,i(a;) 

With HarmonicSums this can be handled as follows. 
HarmonicSums session. 

in[2i]:= RemoveLeadingl[RemoveTrailingO[H[l, -1, 0, x]]] 

Out[2i]= H[0,a;](H[-l,a;]H[l,a:] -H[-l,l,x]) -H[l,x]H[0,-l,x] +H[0,-l,l,x] 



3.2 Multiple Harmonic Sums at Infinity 

Later we will use the fact that the values of harmonic polylogarithms at a; = 1 arc 
related to the values of multiple harmonic sums at infinity. In order to make this more 
precise, we first take a closer look at harmonic sums at infinity. 

Of course, not all multiple harmonic sums are finite at infinity, since for example 
lim^^oo -S*! (n) does not exist. In fact, we have the following lemma, compare [27]: 

Lemma 3.2.1. Let 01,02, •••fip G ^/{O} for p G N. The multiple harmonic sum 
Sai,a'2,...,ap{n) is Convergent, when n — > 00, if and only if ai 7^ 1. 

Similar as for harmonic polylogarithms we say that a harmonic sum 501,02, 
leading ones if there is an i G {1,2, ... ,p} such that for all j G N with 1 < j < i, 
aj = 1. In a similar way as for harmonic polylogarithms we can use (2.1) to extract 
leading ones; for the next example compare [5, Equ. 3.23]. 

Example 3.2.2. For n G N, 

Si,i,2(n) = Si,i{n)S2{n) + ^S^n) + 5i(n)(53(n) - 52,i(n)) - 83,1(71) 
+S'2,i,i(n) + ^54(77.) - S2,2{n) 

Again this can be handled with HarmonicSums: 
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HarmonicSums session. 

in[22]— SRemoveLeadingl[S[l, 1, 2, n]] 

o,t[22i= ^^^'"J"^^^^'''^ + S[l, n] S[3, n] + - S[l, n] S[2, 1, n] - S[3, 1, n] + S[2, 1, 1, n] 



We will always be able to express a multiple harmonic sum as an expression consisting 
only of sums without leading ones and sums of the type Si^. (n) with ttj € N. Since sums 
without leading ones are convergent, they do not produce any problems. For sums of 
the type Si^{n) we have the following proposition which is a direct consequence of 
Corollary 3 of [16] and Proposition 2.1 of [23]. 

Proposition 3.2.3. For w>l and n G N, we have 

Sijn) = j^i^iy-' (fj I = -^YU. . . , (i - l)!S,(n), . . .), 
i=l ^ ' 

where Y-wi^- . . , Xj, . . .) are the Bell polynomials. 

Note that the explicit formulas were found empirically in [11]; for a determinant eval- 
uation formula we refer to [11, 15, 2]. Using this proposition we can express each sum 
of the type Si^{n) by sums of the form Si{n), where z G Z. Hence we can decompose 
each multiple harmonic sum Sai,a2,-,ap{n') in a univariate polynomial in Si{n) with 
coefficients in the convergent harmonic sums. So all divergences can be traced back to 
the basic divergence of <S'i(n) when n ^ oo. 

Notation 3.2.4. We will define the symbol S'i(oo) := lim„^oo "S"! (n) , and every time 
it is present, we are in fact dealing with limit processes. For k & N with k > 1 the 
harmonic sums Sk{oo) turn into zeta-values and we will sometimes write for S'fc(oo). 

Example 3.2.5. 

lim 5i,i,2(n) = lim \ls1infS2in) + -Si(n)(53(n) - S2,i(n)) 
- -53,1(71) + ^SAin) + 52,i,i(n)| 

= 5i(oo)2|-5i(oo)C3 + ^ 

For the computation of the actual values of the harmonic sums at infinity see [41]. 
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3.3 Values at One and Power Series Expansion of Harmonic 
Polylogarithms 



Due to trailing zeros in the index, harmonic polylogarithms in general do not have 
regular Taylor series expansions. To be more precise, the expansion is separated into 
two parts, one in x and one in log(a:). E.g., it can be seen easily that trailing zeros are 
responsible for powers of log(a;) in the expansion since [33]: 



lo,{x) = ^^j\^\og\x) = x-+'Y. 



Subsequently, we only consider the case without trailing zeros in more detail. For 
weight one we get following the well known expansions. 

Lemma 3.3.1. For x G [0 , 1) , 

oo j 

Hi(x) = -log(l-a;) = ^^, 

i=l 

oo 

H_i(x) = log(l + a;) = ^ 



i=l 



For higher weights wc proceed inductively. If we assume that m has no trailing zeros 
and that for x G [0, 1] we have 



1=1 



then for a; G [0 , 1) the following theorem holds. 
Theorem 3.3.2. For x G [0 , 1) , 



oo j 

Hl,m(a;) = ^—Saa,n{'i'-'^)=^—S<Ta,n{'i')-^^^Sn{i), 



1=1 

oo oo V oo V 4 



% ^ — ' I — ' Z" 

i=l j=l 1=1 

oo 



1 ' • 1 ' • 1 ' 

1=1 1=1 2=1 



In order to prove this theorem, we need the following lemma: 
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Lemma 3.3.3. Let a G N,(7 G {—1, 1} and n be a vector with integer entries. Then 

Sn{i + 1) 



lim 

i— >-oo 



and the power series 



Snii) 



1, 



is uniformly convergent for x G (0, 1). 



Proof. If n docs not have a leading one then Hrrii^ooSnii) exists, so it is clear that the 
first statement holds. On the other hand if n has leading ones, we can first extract 
them, and we are only left to show the statement for n= (1) : 



lim 



Si{i + 1) 



lim „ , , ^ ^ = lim 



+ 1 = 1. 



Si{i) i->oo Si{i) i-^oo {1 + i)Si{i) 
In order to prove the second statement we apply the ratio test [18, p. 205]: 

c'+lc7'+^9„(j+l) 



xSn{i + l)i" 



{i + \YSn{i) 



X 



i" 


Snii + 1) 




Sn{i) 







X. 



Hence using the criterion of Weierstrass [18, p. 555] the power series converges uniformly 
for xG (0,1). □ 



Now we can give the proof of Theorem 3.3.2. 



Proof. Prom the definition of harmonic polylogarithms and the theorem of dominated 
convergence of Lehesgue (see e.g. [40]), which we can apply due to Lemma 3.3.3 and 
which will allow us the exchanges of sum and integral, we get: 



Ho,m(^) = / -Y.-^^r^{^)dy 
Jo y * 

i=l ^ 

OO j j 



■ 1 * 

1=1 
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^ i=l "^'^ \fe=0 / \i=l / 

^ OO {*X 

= / Ey'^''^-'^.»(^ + i)^y= / Ey'^^-.-w^?/ 

i=\ i=0 

OO 

= ^ ] ~ Saa,n{i ~ !)• 



a; 

The third part fohows analogously to the second part. □ 
Example 3.3.4. For x G [0, 1] , 

OO OO 



i=l i=l 



With HcLrmonicSums this job can be done as follows. 



HarmonicSums session. 




In[23]:= HToS[H[-l, 1, x]] 




Tj r ^ ^ {-xyS[-l,i] 




ln[24]:= HToS[H[-l, 1, 0, 1, x]] 













If we want to analyse the values of harmonic polylogarithms at x = 1, we can first 

extract trailing zeros. According to Remark 3.1.15 we get a polynomial in log(x) with 
coefficients in the harmonic polylogarithms without trailing zeros. If wc send x to one 
and if the coefficients are continuous for x € [0, 1] then only the constant term (the 
term without a power of log(a;)) will remain, since log(l) = 0. Note that he coefficients 
may not be continuous for x € [0, 1] since they can contain harmonic polylogarithms 
with leading ones. However, according to the following lemma even the terms with 
coefficients which contain harmonic polylogarithms with leading ones will vanish if 
they contain a power of log(x) and if we send x to one. Hence we only have to look at 
those functions without trailing zeros. 
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Lemma 3.3.5. For p, g € N, 

lim logP(x) log'?(l -x)=0 



Proof. The proof follows by de I'Hospital's rule (see e.g. [18, p. 287]) and induction. □ 



As worked out earlier, the expansion of the harmonic polylogarithms without trailing 
zeros is a combination of sums of the form: 

oo 



ae{-l,l}, a EN. 



X 

1=1 

For X ^ 1 these sums turn into harmonic sums at infinity if era 7^ 1: 



00 



i= 

If era = 1, these sums turn into 



1=1 



00 

.1 



E 

Sending x to one gives: 



i=l ' 



lim > x = 00 

x^l ^ I 

i=\ 

Wc sec that these limits do not exist: this corresponds to the infiniteness of the harmonic 
sums with leading ones: limfc^oo >S'i,n(^) = 00; see Lemma 3.2.1. Hence the values of 
harmonic polylogarithms at one are related to the values of the multiple harmonic sums 
at infinity. 

Example 3.3.6. 

H_i,i,o(l) = -253(00) + 5_i,_2(oo) + 5_2,-i(oo) 
Using HarmonicSums we can carry out this evaluation at 1 automatically. 



HarmonicSums session. 

in[25]— RemoveTrailingO[H[-l, 1, 0, x]] 

Out[25i= H[0, x\ H[-l, 1, x\ - H[-l, 0, 1, x\ - H[0, -1, 1, x\ 

ln[26]:= % /. H[0, x] ^ 

Out[26]= -H[-l, 0, 1, a] - H[0, -1, 1, x\ 
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ln[27]:= % // HToS 



/ OO ^ \ oo 
Out[271= -2 (^XJ Jj+E 



(-ccrS[-2,^] 



oo 



(-:c)'S[-l,»] 



We send x to 1. Afterwards we rewrite the output in terms of harmonic sums, this is done by 
using the function TransformToSSums. 

in[28]:= % /. X 1 // TransformToSSums 

Out[28]= -2S[3,oo] + S[-2,-l,oo] + S[-l,-2,oo] 

Wc give a slightly bigger example: 

in[29l:= HToS[H[-l, 0, 1, 0, 0, 1, x]] /. X 1 // TransformToSSums 
Out[29]= -S[6, oo] + S[-l, -5, oo] + S[3, 3, oo] - S[-l, -2, 3, oo] 

Remark 3.3.7. Details on the procedure TransformToSSums can be found in Chap- 
ter 5. 

There is one additional aspect, which causes difficulties when values at x = 1 are 
considered. As mentioned earlier (sec Remark 3.1.10) the product formula for harmonic 
polylogarithms (3.3) does not hold if one of the objects is divergent. So we have to be 
careful when we work with polylogarithms at x = 1 which are divergent. We state the 
same example as in [33] to illustrate this problem (the subleading terms will be wrong) : 
Prom Lemma 3.3.1 we have 




oo 



i=l 



using (3.3) (here we require that x < 1) and Theorem 3.3.2 we get 




Sending x to one we get 



Hi(l) 



lim > — = 51(00), 




2Hi,i(l) 



lim 2^x^ - 1) = 2(5i,i(oo) - ^2(00)) 



i=l 



iSi{oo)f - 82(00). 



We see that 



(lim Hi(x))2 ^ lim(Hi(a;))2. 
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We can however use multiple harmonic sums to solve this problem, since there the 
product also holds for sums at infinity. Namely, we can replace a sum at infinity by 
a sum with upper bound N where N is large but finite. Then we can use the shuffle 
algebra to cancel possible divergencies and afterwards we send N to infinity. Following 
this approach we get 

(lim Hi(x))2 = (Si(oo))^ = 2Si 1(00) - ^2(00) = 2 lim Hi i{x) + lim H2(x). 

This way is consistent and it will allow us to define the McUin transform properly 
in the next section. Summarizing, we can again express all divergencies by the basic 
divergency 51(00), which is equal to lim^j^i Hi(a;). 



3.4 The Mellin- Transform of Harmonic Polylogarithms 



In order to accomplish differentiation of harmonic sums, we look at the so called Mellin- 
transform of harmonic polylogarithms; compare [30]. 

Definition 3.4.1 (Mellin- Transform) . Let f{x) he a locally integrable function on (0, 1) 
and n G M. Then the Mellin transform is defined by: 

M {f{x),n) = / x'^f{x)dx, when the integral converges. 
Jo 

Remark 3.4.2. A locally integrable function on (0, 1) is one that is absolutely integrable 
on all closed subintervals 0/ (0, 1). 



For f{x) = 1/(1 — .x) the Mellin transform is not defined since the integral Jq does 
not converge. Like in [33] we extend the definition of the Mellin transform to include 
functions like 1/(1 — x) as follows 

Definition 3.4.3. Let f{x) be a locally integrable function on (0, 1) and let g{x) be a 
harmonic polylogarithm which is finite at x = 1, and let ^? G N. Then we define 

M+(/(x),n) = M(/(a!),n) = j x^fixjdx, 

\l-x J Jq 1-x 
^+ / g(x)log^(l-x) \ ^ {x-g{x)-g{l))logP{l-x) ^^ 
V '/ Jo 

when the integrals converge. The function f is supposed to be a harmonic polylogarithm 
without leading ones when the factor 1/(1 — x) is present. 



Harmonic Polylogarithms and the Differentiation of Harmonic Sums 



44 



According to the definition, if we want to compute the Mellin-transform of liarmonic 
polylogarithms, we have to take care of harmonic polylogarithms which are not finite 
at X = 1. As mentioned in Remark 3.1.1 only harmonic polylogarithms with a leading 
one (and not followed only by zeros) are not finite at x = 1. By using equation 
(3.4) we can extract all the powers of log(l — x) which cause the infiniteness at a; = 
1. After this extraction the remaining harmonic polylogarithms are finite at x = 1. 
The following lemma will guarantee the existence of the Mellin transform of harmonic 
polylogarithms. 

Lemma 3.4.4. For a harmonic polylogarithm Hj„(x) the integrals 

/ X iim{x)dx and / — — 

Jo Jq 1 + x 

converge for x G [0, 1] and n G M. //Hj„(x) does not have leading ones and p G N then 

x^H^(x)-H^(l) ^^^^^ (x"H^(x) - H^(l)) log^(l - x) ^^ 
Jo ^-x Jo 1 - a; 

converge for x G [0, 1] and n G M. 



Proof. Let n G M and consider x"Hm(x)(ix. After extracting trailing zeros and lead- 
ing ones we get a combination of integrals of the form: 

/ x" log"(x)Hyfi(x) log''(l - x)dx, a,be No, 
Jo 

where Hm(.T) is continuous on [0,1]. It suffices to prove, that these integrals converge. 
|Hrn(x)| is continuous on [0, 1] and hence has a maximum C on [0, 1]. Thus 

[\''log%x)Rru{x)log\l-x)dx < [\''\\og%x)\\Rrn{x)\ l0g\l-X 

Jo Jo ^ V ' 



dx 



<c 



< 



C [ x'*|log"(x)| log''(l-x) 
Jo 

1 

limC x" |log"(x)| log''(l-x) 

•^f <1 



dx 



+ limC 

e->-0 



<|logi|''=:Ci 



X 



<1 



, |log'^(x)| log^(l-x; 

V ' 

<|logir=:C2 



dx 



< CCilim / |log"(x)|dx 



+CC2lim / 



log''(l - x] 



dx. 
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Since the last two integrals eonverge, we finished this part of the proof. 
The convergence of the second integral Jq ^ ^mi^)dx proved analogously. From 

now on we assume that Hm(x) does not have leading ones, hence H^(l) is finite. If we 
find xi G M with < xi < 1 such that for all x G (xi, 1) we have 



x'^H„(x) - H„(l) 



1 -X 



<{l-x)-2, (3.5) 



the convergence of /o is proved; this follows since the integral 

(1 — x)~^ converges. Since x"Hj„(x) — 11^(1) is continuous on [0,1], there is X2 
with < X2 < 1 such that x"''H.m{x) — 11^(1) is either < or > on [x2, 1]. We will 
prove the > case, the other case follows analogously. For x G (x2, 1) we want to show 
(3.5) or equivalently: 

x'^H^(x) - H„(l) 



X) 2 



< 1. 



Let m = (mi, 1712, ■ ■ •) with mi G {0, —1} by assumption. With de I'Hospital's rule, see 
e.g. [18, p. 287], we get: 

X"H„(X) - H^(l) ^. nX^-^Rmix) - j^nm,,ms,..Xx) 

lim — J — — = lim ^ 

^-^1 (l-x)2 _i(i_ar)-2 

= lim (-2nx^-^{l - x)^H^(x)) 



+ lini ( 2—^ (1 - x)2H^2,m3,...(2:) 



Since Hm(l) is finite, the first limit exists and is zero. If 7712 7^ 1 then Hm,2,m3,...(l) exists 
and so the second limit is finite and equals zero. If 777-2 = 1, we can extract leading ones 
and split the limit. We will get limits of the form 

lim U^^il - x)i log (1 - x)'=H„,...(x)^ 
x^i \ X — mi J 

where Ha^...(x) is finite at 1 and A; G N. Since lima;_;.i(l — x)^ log (1 — x)^ = 0, all these 
limits vanish. Hence there is xi G [x2, 1] such that (3.5) holds and so '^"H"-(g-H„.(i) 
converges. 



For the last integral we want to find again xi such that 

(x"H^(x)-H^(l))log^(l-x) 

(1 — X)2 

holds for all x with xi < x < 1. Again we can find X2, with < X2 < 1 such that 
(x"H„(x) - H^(l)) log^'(l - x) is either < or > on (x2, 1). 
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We will prove the > case, the other case follows analogously. Using again de 
I'Hospital's rule we get 

lim (^"^"'(^^ ~ ^"g''^^ ~ ^) = 

^-^1 (1-,t)5 

hm (2(1 x)- ^(^"""-(^) - ^-(^)) ^°g'"'(^ - 

x^l V 1 - X 

- lim (2(1 - x)^ log*^ (1 - x) fnx"-iH^(x) - -^^H^,,„3,...(x)') ) . 

For the second limit we can use the same arguments as we used in the proof of the 
convergence of the third integral to show that its limit is zero. For the first limit we 
can use induction on p to see that the limit is zero too. Hence there is xi G [x2, 1] such 
that (3.6) holds and so (^"H"'(a;)-Hm(i)) iogP(i-a;) converges. □ 

Remark 3.4.5. Prom now on we will call the extended Mellin transform M+ just 
Mellin transform and we will write M instead of M'^. 

Example 3.4.6. For x G [0, 1] and ra G M 



M(Hi,o(.x),n) = y"'x"Hi,o(x)dx, 
M(^,n) = j^^"Ho,(x^-Ho,(l)^^^ 

Mi^,n] = ^ f Ho(x)H,(x) _ HcuM 



1 — X / \ 1 — X 1 — X 



1 — X 'y \1 — x' 

(x'^Ho,i(x)^Ho,i(l))log(l-x) ^^ _ x"Ho,i(x)-Ho,i(l) ^^_ 
Jo 1 - a; Jo 1 - ^ 



compare also [11]. 



Now we want to study how we can actually calculate the Mellin transform of harmonic 
polylogarithms. We will proceed recursively. It is quite easy to get the Mellin transforms 
of harmonic polylogarithms with weight one [33]; see also [11]. 

Lemma 3.4.7. For n G N we have 

M(Ho(x),n) = ^ 
M(Hi(x),n) = 



(n + l)2' 
Si{n + l) 



n + 1 ' 

M(H_,w,„) = (-ir^=;f^ + 1f^(i + (-i) 
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Proof. We just give a proof of the first identity. By using partial integration we get for 



/•I n+1 ^ /•! n+1 

M{Roix),n) = / x''Ro{x)dx = ^Ro{x) -/ ^ 



(n + l)2 



(n + l)2 



□ 



The higher weight results can now be obtained by recursion; see [33, 11]. 
Lemma 3.4.8. For n G N and m G {0, -1, 1}^, 

77. + 1 n + 1 Jo 



M(Hi,^(a;),n) = 



Hl,m(l) , 1 



+ 



n + 1 n + 



{x)dx I , 



M(H_i,^(x),n) = H-i^(l) a;"H^(x)da; - -^H^(x)da; ) . 



Proof. We get the following results by partial integration: 



/ x"Ho,m(x)da 

JO 

/ x''Hi,^(x)da 



JO 



„n+l 



-Ho,m(x) 

n + 1 

Ho,m(l) 1 



1 



^0 

1 



n + l 



H.m{x)dx 



n + l n + l Jq 



x"'H^(x)dx, 



lim 



X = lim / x"Hi^^(x)(ix 



„n+l 



e^l V 77. + 1 



Hi,m(a;) 



„n+l 



io (n + l)(l-a;; 



-Hm(x)dx 



^H_,„(x) ^ - ^-^^^^H„(x)dx 

^^^^ - ^ ^'^(1 - :^)H^(x)c/x. 
n + l n + lJo l + x 



□ 



Remark 3.4.9. We want to point out again that Hi^m(l) '■= lima;-).i H\^m{x) and hence 
we are dealing with limit processes. 
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As wc could sec in the previous lemma, we will need the results for harmonic polylog- 
arithms weighted by 1/(1 + x) or 1/(1 — x). We start with harmonic polylogarithms 
with weight one; see [33, 11]. 

Lemma 3.4.10. For n G N and x G (0, 1), 

'Ho(x) 



M 

M 

M 

M 
M 
M 



1-x' 

Ho(x) 
1 + x 

H-i(x) 
1 — x 

H-i(x) 
1 + x ' 
Hi(x) 
1-x' 
Hi(x) 
1 + a; 



,n 



n 



n 



,n 



= 52(n)-C2, 
= (-ir(-|-5_2(n)), 



log'(2)-C2 



2 

n / V(2) 



= (-1) 



= --51,1(71), 



+ log(2)(5_i(n) - Si{n)) + 5-i,-i(n), 
+ Iog(2)(5_i(n) - Si{n)) - -Si,_i(n)^ , 



(-If 5_i,i(n) + 



C2-log'(2) 



Proof. We use again partial integration and Lemma 3.4.7. In addition we use tables by 
Vermaseren (see [41] and [33]) where we can find values of multiple harmonic sums at 
infinity and harmonic polylogarithms at one. To verify the exchange of integrals and 
summations we used the theorem of dominated convergence of Lebesgue (see e.g. [40]): 

\l-x J Jo 1-x Jo 

00 „1 00 J 

= 5Z / ^'Ho(x)da; = - ^ 

i=n ^ i=n 
00 J 

= - ^ ^2 = S2{n) - 82(00), 

00 „i 00 /_iy+l 

= (-irE/ (-^rHo(x)dx = (-i)'^j:i-^ 



E ^ = (-ir(5-2(oo)-5_2(n)), 

j=n+l 



1 — X ^ ) Jq 1 — X 

,1 00 



/•I ^ 

/ 

•''0 tZ 



x*H_i(x) -a;*-"H_i(l))dx 
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oo „i 

V/ (a:^H_i(ar)-x^-"H_i(l))da; 

oo ^ „1 „i ^ 

x'B.-i{x)dx - j a;*-"H_i(l)da; 

p-i(i + l)(-ir + H_i(l)((-l)- + l) _ H_i(l) \ 
N^oo -^V i + 1 i- n + l) 

i=n ^ 

I Y ^^^(^ + + H-i(l)((-ir + 1) H_i(l) \ 

\i=n 1=1 / 

lim (5_i,_i(n) + H_i(l)(5_iH - - 5_i,_i(iV + 1) 

-H_i(l)(S-i(Ar + 1) - Si{N + 1)) -H_i(l)5i(Ar - n)) 
5_i,_i(n) + H_i(l)(5_i(n) - Si(n)) - -S_i,_i(oo) - H_i(l)S_i(oo) 

^"^'^^j"*^' + log(2)(5_i(n) - Si(n)) + 5_i,_i(n), 



(1-x)' 



oo »i 



i=n 
oo 



, .nY- .g-i(i + l) + H_i(l)((-l)- + l) 

= i-ir lim (5i,_i(iV) + H_i(l)(5i(iV)-5_2(A^))-5i,_i(n) 
-H_i(l)(5i(n)-5_2(n))) 

= (-ir +log(2)(S_iH - 5i(n)) - , 



= lim > 

V-j-oo-^ 

i=n 

lim (^i 

f— >(X) 



Ar->-oo ^-^ \ i + 1 i — n + l 

i=n 

lim (5i,i(iV + 1) - Si,i(n) - 5i,i(iV - n + l)) 



( i).f- (-i)''^'(' + i) . ( 1). f; (-^)'Si(i) 

i=n i=n+l 

= -(-l)"(5_i,i(oo)-5_i,i(n)) 

. (-ir(s_.,(„) + «^^). 



□ 
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The higher weight results can be again obtained by recursion, however we will first 
show uniform convergence for a series which will pop up. 

Remark 3.4.11. In the following recursions we will use the following definitions: For 
n G N and x G (0, 1), 

5Q(n) := S{n):=l, 
H()(x) := H(x):=l. 

Lemma 3.4.12. Let m be a vector with components in {0, 1, —1}, p be the empty vector 
or a vector with components in the non zero integers and let k E NU {0} . Then the 
series ^ 

converges uniformly for x G (0, 1). 



Proof. We apply the ratio test [18, p. 205] and use Lemma 3.3.3: 



lim 

i— >-oo 



X nm{X) (j^_2)fc 



X rirriyx) (^^.^i^k 



lim 



= X lim 

i—^oo 



Spii + 2) (i + 1)* 



Sp{i + 1) (i + 2y 



Sp{i + 2) 



Sp{i + 1) 



lim 



(Z + 1)* 



(i + 2)* 



x. 



Hence using the criterion of Weierstrass [18, p. 555] the series converges uniformly for 
a;G(0,l). □ 



In the following three lemmas we will present the recursion (compare [33]), which we 
will use to compute the Mellin-transform of general harmonic polylogarithms. 

Lemma 3.4.13. For x G (0, 1), n G N, /c G N U {0}, m a vector with components in 
{0, 1, —1} or m = and p a vector with entries inN or p = {) we have: 

ri °^ S* (i + 1) 

J ^i-^y^o,m{x)-^-y^^dx = -Ho,m(l)('S'-A;-i,p(oo) - 



(i + l)'' 

■■'^ Spii + 1 

L) 



i=n 
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Q (j -L -[) 

/ ^ (-x)'H_i,^(x) '.yj dx = -H_i,„(l)(5_fe_i,p(oo)-S_,_i,p(n)) 

„1 oo 

-/ Vx^H4x)(-(-ir5fe+i,p(z+i) 

•^0 i=n ^ 



Proof. We give a proof of the first identity. The second follows similarly. By Lemma 
3.4.8: 

^ i=n ^ ' i=n ^ ^ ^ 



1 ■ 

- T—r / x*Hm(x)(ia 



-Ho,m(l)('S'_fe-i,p(oo) - S-k-i,p{n)) 



□ 



Lemma 3.4.14. For n & N, k £ NL) {0} , m a vector with components in {0,1, —1} or 
m = and p a vector with entries inN or p = {) we have in the ring of formal power 
series: 

/•I °° C (A _|_ I) 

I ^{-xyR,,m{x) rr^J dx = -Hi,„(l)(5_fe_i,p(oo)-S_,_i,p(n)) 

^ i=n ^ ^ 

1 oo 

- / Vx'H„(x) -S_fe_i,p(z + 1) 

•^0 i=n ^ 



Lemma 3.4.15. For n eN, k e NU {0}, m a vector with components in {0, 1, —1} or 

m = {) and p a vector with entries inN or p = {) we have in the ring of formal power 
series: 

/ ^^'^0,rn{^)-fr—V^dx = Ho,m(l) ('S'a;+1,p(oo) - S'fe+i,p(n)) 



-n 



dx, 
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1 °° . S + 



'0 i=n 



/•l oo 

- / y2x'Hmix)[Sk+l,p{i + l) 

i=n ^ 

Spji + 1) ^ , A , 



1 oo 



(1 + Ij"*"^^ ' ' 

(3.7) 



Proof. We only show the third identity. It is imphed by 



oo 



+ V ^ + 1 Jo ^ + 1 



= H_i,m(l)(5'A;+i,p(oo) - S'fc+i,p(n)) 
= H_i,m(l)(5fc+i,p(oo) - Sk+i,p{n)) 

1 oo 

-/ Vx*H^(x)((-ir5_fe-i,p(i + l) 
+^^.-(-'yS->^-^An))dx; 

in the last equation we used: 

i=0 ^ ^ j=Q i=0 j=0 ' 



£(-i)VH„(x)i:(-iy^|yjtl) 

i=o j=o ^-^ 
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i+l 



^(-l)VH^(x)^-(-l 



i=Q 



k+1 



^ (-l)VH„(x)5_fe_i,p(i + 1). 



i=0 



□ 



Remark 3.4.16. In Lemma 3.4-13 the left sides and the right sides converge and 
hence the equalities hold analytically. In Lemma 3.4-14 the left side converges, while 
the convergence of the right hand side is not obvious since Hi^m(l) is present. In Lemma 
3.4.15 the left sides do not converge and so do the right hand sides. 
The left hand sides of 3.4.13 - 3.4.15 can appear in the compnitation of the Mellin- 
transform and hence we use the equalities formally. It turns out that the divergencies 
cancel in the end and all divergent parts disappear. At this point one would have to 
prove that the Mellin-trans forms and the results from these recursions are equal in an 
analytic point of view, however this is omitted here. 



We are now able to compute Mellin transforms of all harmonic polylogarithms. We can 
use the previous lemma even if neither 1/(1 + x) nor 1/(1 — x) is present, since 



1 + a; 1 + X 

Hence we do not need Lemma 3.4.8 to compute the Mellin transform of harmonic 
polylogarithms (but we used it to establish the Lemmas 3.4.13 - 3.4.15). In the base 
case one has to evaluate harmonic polylogarithms at x = 1 and harmonic sums at 
infinity; for this task we use tables from [33] which contain expressions up to weight 9. 
The results of the Mellin transforms should be finite, hence the divergencies introduced 
in the Lemmas 3.4.13-3.4.15 should cancel. The result of the Mellin transform is an 
expression with multiple harmonic sums in the argument n. The algorithm can be 
summarized in Algorithm 3.4. We illustrate the usage of the previous lemmas on two 
examples: 

Example 3.4.17. For n eN, we perform the following formal calculus. 

^(iT?''') = (-irj['f;(-^rHi(x)dx = (-ir^'f;(-xm(x)|^dx 

„i 00 

= -(-irHi(l)(S_i(oo)-5_i(n))-(-ir / J2^\-S^,ii + 1) 

^ i=n 
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{i + iy 

I oo 



-(-l)'^Hi(l)(5_i(oo) - 5_i(n)) - (-1)'^ / V -x'S-iii + l)da 

•^0 i=n 



(-1)" (-Hi(l)(5_i(oo) - S_i(n)) + 5i,_i(oo) - Si,_i(n) - 5_2(oo) 
+5_2(n) + S-i(n)(5i(oo) - 5i(n))) , 

1 - a; y ./n 1 - X ./n 



„1 oo 

-H_i,o,i(l))dx 

„1 / oo oo \ 

/ ^x^H_i,o,i(x) - H_i,o,i(l)i;xMdx 

\i=n i=0 / 

„1 oo „i oo 

/ Va;^H_i,o,i(x)da;- / H_i,o,i(l) V x'dx 

„1 oo 

/ ^a;*H_i,o,i(x)dx - H_i,o,i(l)Hi(x) 

^ i=n 

C2-log'(2)~ 



(-ir 5_i,i(n) + 



see Lemma 3.4-10. 



Finally let us look at the result for a bigger example: 



Example 3.4.18. 

+ nuJ^) I I 3C2log(2)^5-i(n) 



2y 20 4 



6 V2. 
33C|gi(n) 3C2log(2)^>gi(n 
20 4 



Harmonic Polylogarithms and the Differentiation of Harmonic Sums 



55 



Algorithm 1 Mellin Transform 



input: fix) = Rmix), f{x) = ^ or f{x) = %f 
procedure MEhhlN {f{x),n) 
r = 

if f{x) = H„(x) then 

s = Mellin(^ 
else if f{x) 



'IM. 

l+x 



nj + MELLlN({^,n + 1) 
then 



else if f{x) = and m[l] ^ 1 then 

else if f{x) - 

'^^ \ T*" I— I I /-y I 



and m[l] 



1 then 



X Hm(x) c,;^_^-|^\o'' dx 



* Jo ^i=n 

r = REMOVELEADINGl(Hm 

nsum =number of summands in r 
modify r in the following way: 
for i = 1 to nsum do 

if the ith summand of r equals c * Hi^i^...^i(x) then 

add 1 to the index set and set x = 1 
else if it equals c * Hy(x)Hi^i^...^i(a;), y 7^ (1, 1, . . . , 1) then 

add 1 to the index of Hij....^i(x) and set x = 1 
else if it equals c * }ly{x), y / (1,1,...,!) then 
multiply Hi(l) to the summand and set x = 1 
end if 
end for 
end if 



> Use 3.4.13-3.4.15 

> first index 7^ 1 

> Use 3.4.13-3.4.15 

> Use 3.4.13-3.4.15 

> first index = 1 

> Use 3.4.13-3.4.15 



> c G 



s = s — r 

transform all Hy(l) in s to Sy{oo). 
expand all products in s 
return s 
end procedure 



Use Theorem 3.3.2 
> Use (2.1) 
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+ log(2)5_,_,_.(n) + ^%^ 

- log(2) S_i,_i,2(n) + 5_i,_i,_2,-i(n) + ^ 
13 1og(2)2C3 13 1og(2)5_i(n)C3 



+ 



8 4 
13 1og(2)5i(n)C3 845 Cs 
4 64 



The Mellin transform of harmonic polylogarithms is implemented in HarmonicSums: 



HarmonicSums session. 

in[30):= Mellin [H[-l, 1, 0, x]/(l + x), x, n] 

Out[30]= (-1)" (H[-l, -1, 1, 0, 1] + H[-l, 1, 0, 1] S[-l, n] - H[-l, 1, 0, 1] S[l, n] + S[l, -1, 2, n]) 

in[3i]:= % // ReplaceByKnownFunctions 

/ ,^ ln2'' In22z2 722^ 131n2z3 ln2z2S[-l,n] , , ln2z2S[l,n] 

Outi3ii= (-1)" (21i4half+^^ — + + _L^_z3S[-l,n] 7T^^ + 



z3 



S[l,n]+S[l,-l,2,n]) 



in[32]:= Mellin[H[-l, 1, 0, -1, x]/(l - x), x, n] // ReplaceByKnownEYinctions // Expand 

■,-,,-r, ,r o,-., ,r, o ln2'"' 63 ln2 z22 131n22z3 3z2z3 845z5 , , 

Out[32]= Illi5half+31i4halfln2+ H — 1 — 41i4half S[-l, n] - 

ln2^S[-l,n] 31n2^ z2S[-l,n] 33z2^S[-l,n] _ 13 ln2 z3 Sf-l, n] 

6+4+ 20 4 + 

ln2''S[l,n] 31n22z2S[l,n] 33z22S[l,n] 131n2z3S[l,n] 

Q 4 w + 4 + ln2b[-l,-l,-2,nj + 

z2S[-l,-l,l,n] _ s[_i^ 2, n] + S[-l, -1, -2, -1, n] 



Remark 3.4.19. The procedure ReplaceByKnownFunctions uses tables by Vermaseren 

to replace harmonic sums at infinity and harmonic polylogarithms in one by known 
functions. For these known functions we use the notation used by Vermaseren. The 
tables can be found at www. nikhef.nl/^ form. 



3.5 The Inverse Mellin Transform of Multiple Harmonic 
Sums 



We already saw that the Mellin transform of harmonic polylogarithms can be expressed 
in terms of multiple harmonic sums. Now we want to go the other way round. 
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In order to get the inverse Mellin transform of multiple harmonic sums, we have to con- 
sider distributions. In fact we will consider distributions which are either differentiable 
functions in the class of C°°((0, 1)) or are Dirac-(5-distributions 6{1 — x) G D'[0, 1]. We 
want to point out that the inverse Mellin transform of constants is 5{x — 1) where 6{x) 
is the Dirac-(5-distribution and that d{x — 1) will only appear as the inverse Mellin 
transform of constants. Since we are mainly interested in the differentiation of multiple 
harmonics sums and since a constant vanishes if it is differentiated, we will skip further 
details on the inverse Mellin transform of constants. Subsequently, we will write M~^{c) 
for the inverse Mellin transform of a constant c. 

To prepare the stage we give some properties of the shuffle product of harmonic sums; 
here the product among harmonic sums is defined as for harmonic polylogarithms in 
Definition 3.1.6. 



3.5.1 Properties of the Shuffle Product 



Definition 3.5.1. For ak G ^/{O} with 1 < k < I and c G ^/{O}, we define 

^ai,a2,...,ak-i,{ak,c),ak+i,...,ai ('^) 'S'aj^,a2,...,Ofc_i,ajfc,c,a/;+i,...,a; ('T') 

'^Sai,a2,...,ak-i,ak,0'k+i,c,ak+2,---,ai (^) 
+ ■ ■ ■ + Sai,a2,...,ai_2:ai^i,c,ai{n') 

+5'ai,a2,...,ai_i,a;,c('T'); 

similarly, we define 

(ai, 02, . . . , afe_i, {ak, c) , ak+i, ■ ■ ■ ,ai) := {(ai, 02, . . . , afe-i, a^, c, ak+i, ■ ■ ., ai), 

(ai, 02, . . . , ak-i,ak, a^+i, c, ak+2, ••-,«;) 
, • • • , (oi, 02, . . . , ai-2, ai-i,c, ai), 
{ai,a2,...,ai^i,ai,c)}. 



Example 3.5.2. 

5'l,(2,3>,4,5H = -51,2,3,4,5 W + 'S'i,2,4,3,5(n) + S'l,2,4,5,3W- 



Lemma 3.5.3. Let ai, bj G Z/{0}, for 1 < i < k and 1 < j < / and consider the shuffle 
product S'ai,a2,...,afe(^)'-'-''S'6i,62,..-,6i(")- ^ representation of 



'S'ai,a2,...,afc(n)LuS'6i,62,...,6,(n) 
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such that it is linear in the harmonic sums . If we replace every occurrence of bi in 
P by {bi,bi+i) and call the result P, then we get (using Definition 3.5.1) that P is a 
representation of 

We give an example which illustrates this lemma. 
Example 3.5.4. Consider 

,02,61 'S'ai ,61 ,02 ,01,02 

(n). 

If we replace bi by (61,62) we get: 

-^ai, 02,(61 ,62} (") + '5'ai, {61,62}, 02 H + 'S'{6i,62>,oi,02H = S ^a-i M ,b2 i^) + -5*01,61,62,02 ("-) 

+-5*01,61,02,62 (^) + -561,62,01,02 (^) 

+ -561,01,62,02 (^) + -561,01,02,62 (^) 

= -501,02 ('^)i-U'56i, 62- 

Remark 3.5.5. // we define aU b for two vectors a = (ai, 02, . . . , a^) and b = 

(61, 62, . . . , 6/) as the set of all mergings of a and b in which the relative orders of 
the ai and bj are preserved (for example (oi, 02, 61, 62, ... , 6;, as, 04, ... , a^) G aU 6 but 
(a2, ai, 61, 62, . . . , 6;, 03, 04, . . . , afe) ^ a U 6) then the shuffle product is nothing else but 

Sa{n)inSb{n) = ^ Sc{n). 

c6aU6 



We give a proof of Lemma 3.5.3: 

Proof. Let a = (ai, 02, ... , a^) and b = (61, 62, . . . , 6;), let A be the set a U 6, however 
replace each bi in A by (6;, 6;+i) and let B = (ai, 02, . . . , a^) U (61, 62, ... , 6;, 6;+i). Let 
A be the union of the sets which arise after using Definition 3.5.1 for all elements in A. 
It remains to show that A = B. 

Wc start by showing A C B : Take an element x£ A. The relative order of the a.j and 
bj in aU 6 is preserved and the exchange of 6; by (6i,6;+i) does not change the order 
of tti and bj. Due to Definition 3.5.1 6; is left to bi^i for all harmonic sums in A. Hence 
X has to be in B. 

Now we show B C A : Obviously it is always possible to find for a given 
X e B axi element y e A such that after using Definition 3.5.1 x occurs 

^Note that this is always possible by using (2.19) 
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in y (for example let x = (oi, 61, 62, . . . , • • • > Qfe, then take y = 

(ai, 61, 62, . . . , bi-i, {bi, , 02, . . . , ak))- Hence x G A. This completes the proof. □ 

3.5.2 The Most Complicated Harmonic Sum 

In order to get the inverse Mellin transform of multiple harmonic sums, we exploit the 
following order on the set of harmonic sums. 

Definition 3.5.6 (Order on harmonic sums). Let Smi{n) and Sjn^{n) he harmonic 
sums with weights wi, W2 and depths di and ^2, respectively. Then 



For a set of harmonic sums we call a harmonic sum most complicated if it is a largest 
function with respect to -<. 

Example 3.5.7. We have 



If we consider the set of these 4 sums then >S'2,i,i(n) is most complicated. Note that -< 
is not a linear ordering, e.g., 

Si,2,3{n) -A 52,i,3(ri) and 52,i,3(n) 7^ Si,2,3(n)- 
If we consider the set of these 2 sums then -S'i^2,3(?^) and -S'2,i,3(ra) are most complicated. 

In the following we show that in the Mellin transform of a harmonic polylogarithm 
there is only one most complicated harmonic sum. 

Proposition 3.5.8. In the Mellin transform of a harmonic polylogarithm Hm(a;) 
weighted by 1/(1 — x) or 1/(1 + x) there is only one most complicated harmonic sum 
Srnin),i.e., 



where all harmonic sums in t occur linearly and are smaller then Srn{n) . 

Let us consider the Mellin transform of the polylogarithm Hmim(a;) weighted by 
oi' 1^ and let w be the weight of H^^rnix). By looking at the recursions in the Lem- 
mas 3.4.13, 3.4.14 and 3.4.15 we can see that the multiple harmonic sums appearing in 
the Mellin transform have weight less or equal w + \. We consider now those harmonic 



S'mi(n) -< Sm2{n), if Wi<W2, 

Smi{n) -< Sm2{n), if wi = W2 and di < ^2. 



Si,i(n) ~< Si,2{n) ~< 5i,i,i(n) ~< S2,i,i{n). 
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sums in the Mcllin transform which have weight w + 1 and which have the maximal 
depth and so are most complicated. 

Where do these sums emerge from? Again looking at the recursions in the Lem- 
mas 3.4.13, 3.4.14 and 3.4.15 we see that if nii = 1 then we will have to do either 
(depending if there is a factor (—1)* present) the following 3 integrals 



in the next step. Likewise, we have to do three similar integrals if nii = —1. In this 
situation (mi = 1 or mi = — 1) the most complicated harmonic sums can only emerge 
from the first or the third integral; there the weight of the polylogarithm is reduced 
by one and the depth of the harmonic sum is raised by one; in the second integral 
however the weight of the harmonic polylogarithm is reduced by one, but the depth 
of the harmonic sum is not raised. In the situation mi = there is just one integral 
and the depth of the harmonic sum is not changed. So these sums emerge from those 
recursion branches where in each step either a first or a third integral is used, whenever 
the entry in m is 1 or —1. 

How do these sums look like? Let us start with an example. 

Example 3.5.9. Consider the Mellin transform of MmamM^ jj take the branch 
where always the first integral is used, we get the sum 5*1^2,3 (^t-)- ff we take the branch 
where always the third integral is used, we get Si(n) * 5*2 (ra) * S3{n). This product will 
produce sums of highest depth: 






or these 3 integrals 






Si{n) * S2{n) * Ssin) = 5i,2,3(n) + -Si,3,2(n) + S2,i,3(n) + S2,3,i(n) 

+53,i,2(n) + S3,2,i(n) + t. 
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Here t is a linear expression in sums of weight 6 with depth less then 3. 

Notation 3.5.10. Let a = (ai, 02, • • • , afe) and h = (61, 62, . . . , 6;). We denote the 
concatenation product by a - b. Hence a - b = (ai, 02, . . . ,ak,bi,b2, ■ ■ ■ ,bi). 

Definition 3.5.11. Let a = (ai, 02, . . . , a^) be a vector of length k, G Z\ {0} and let 

d,neN. We define 

rev(a) := 

and 

part(a) := {Sb^ (n) • Sb^in) ■ ■ ■ Sb^in) \ {bi ■ b2 ■ ■ ■ ■ ■ ba) = a} . 

Let X = S'bi(n) • Sb2in) • • • 'S'b^(n) be an element o/part(a). We call d the length of x. 

Example 3.5.12. For a= (1,2,3), 

part(a) = {Si,2,3(n), S^,2{n)S^{n), Si{n)S2,^{n), Si{n)S2{n)S^{n)} . 

In order to get all most complicated harmonic sums, we first look at the branch where 
we only use the first integral and look at the index set a of this multiple harmonic sum 
(let a be the sign of this sum). All most complicated harmonic sums have weight and 
depth equal to the weight and depth of Sa{n). Each element in part (a) corresponds 
to exactly one branch where we use in each step either the first or the third integral. 
The sign of each element depends on a and the length of the element. If the length of 
the element is d then the sign is —a (— l)*^. Combining all possible branches we get the 
following sum 

^_l^Length(3;)^_ 

a;6part(o) 

If we work out the sum and linearize the harmonic sums, a lot of cancellation will take 
place, in fact we have the following lemma. 

Lemma 3.5.13. Let a = (ai, 02, ... , ak) be a vector of length k, ai G Z\ {0}, then 

(-l)Length(a;)^ ^ (-l)'=rev(a) + t, 

a;6pajt(o) 

where t is a linear combination of sums of depth less than k. 

Remark 3.5.14. If we know the set 

part(a) = {S'6i(n) • Sb^in) • • • ^^^(n) | (61 • 62 • • ■ ■ • bo) = a} 
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for a = (ai, . . . , o^), we can get the set part(a) for a = (ai, . . . , o^, a/j+i); 

part(a) = {S'5j(n) • S'b2(n) • • • 5b^(n) • S'aj^^i(n) | (61 • 62 • • • ■ • &d) = «} 

^ V ' 

part]^(a):= 

U {Sb, in) -S^in)--- S,^,a,^^ (n) I (61 • &2 • . . . • fed) = a} . (3.8) 
^ ' 

part2(a):= 



Now we can give the proof of Lemma 3.5.13. 



Proof. We proceed by induction on the length A; of a = (ai, 02, . . . , a^). If A; = 1, then 
part(a) = {Sa^in)} and 

2;£part(a) "^-"v-^ 

ai 

Thus Lemma 3.5.13 is true for A; = 1. Let us assume that the lemma holds for k and 
let a = (ai, a2, . . . , a^+i). With the notation given in Remark 3.5.14, 



(•_2)Length(a;)^ _ ^ ^_-|^^Length(3;)^ _|_ ^ ^_-|^-jLength(a;)^ 



a;£part(a) a;£part]^(a) 3;gpart2(a) 



Let us first look at A: 



a;6part(o) 



According to the induction hypothesis we get 

E {-l)^^''stHx)^ ^ (-l)'=rev(a) + 1 

a;6part(o) 

where t is a linear combination of sums of depth less than k. Therefore the sums of 
highest depth in A are emerging from Sa^:_^_^{n)Sa^:,ak-i,...,ai{n). By using (2.1) we get 

A = — (— l)*(5'aj,_^j,afc,...,ai(»^) + Sak,ak+i,:.,ai{n) H h Sai^,...,ai,ak+i{''^)) + *A 

where t^ is a linear combination of sums of depth less than A; + 1. 
Let us now look at B. The sums of highest depth in B have depth A; + 1 and result from 
the first and second sum in (2.1) since in the third sum the depth is reduced. Hence 
we can neglect the third sum and just look at the first 2 sums. This implies that it 
suffices to look at the shuffie product and not at the quasi-shuffle product to detect 
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the sums of highest depth. By the induction assumption, the only harmonic sum in 
Z^a;epart(a) {- 1)^"'^^^^^'^^ X of highest depth is Sa,^,ak-i,...,ai{n). Using Lemma 3.5.3 and 
(3.8) we get that the sums of highest depth in B originate from S'^afc,aA:+i>,aj:_i,...,ai ('^)- 
We obtain 

~ ("-*-) 'S'(afc,ofc+i),Ofc_i,...,ai(^) = ('S'afc,afc+i,...,ai('^) + 'S'afc,afc-i,afc+i,afc-2-,ai('^) 

H h S'afc,...,ai,afc+i(n)) +is 

where is a Unear combination of sums of depth less than k + 1. Summarizing we get 

J2 (_l)Length(.)^ ^ A + B = -{-!)'' Sa,^„a„...,aAn) + tA + tB 
a;£part(o) 

This completes the proof of Lemma 3.5.13. □ 
Example 3.5.15. For a= (1,2,3), 

J2 (_i)Length(x)^ ^ -Si,2,3(n) + Si,2(n)53(n) + 5i(n)52,3(n)-5i(n)52(n)S3(n) 

a; 6 part (o) 

= -'S'3,2,i(n) + 5'5,i(ra) + 53,3(71) - Se{n). 

Using the considerations above we get that Proposition 3.5.8 holds with Smin) = rev(a). 
The harmonic sum Sjnin) originates from the branch of the recursion where always the 
third integral is used. 

We emphasize once more that we are mainly interested in the computation of a harmonic 
polylogarithm for which a given harmonic sum is the most complicated harmonic sum 
in the Mellin transform. But first we go the other way round. We have to determine 
rev(a). 

3.5.2.1 Finding the Most Complicated Harmonic Sum 
Given: a harmonic polylogarithm Ilm{n). 

Find: the most complicated harmonic sum Sa{n) and its sign cr in M ^ ^^^^ , ■ 

Of course, we could just compute the Mellin transform and look for the most compli- 
cated sum in the result. However, in our considerations it is crucial that this job can 
be attacked directly. In general, this is possible by applying the Lemmas 3.4.13, 3.4.14 
and 3.4.15 recursively. The most complicated harmonic sum originates from the branch 
of the recursion where always the third integral is used. Looking at the recursions in 
the Lemmas 3.4.13, 3.4.14 and 3.4.15 we see that there are basically two parts of the 
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recursions: the first does not contain the factor (—1)*, and the second contains this 
factor. If the harmonic polylogarithm is weighted by 1/(1 — x), then we are in the first 
part, otherwise if there is the factor 1/(1 + x), we are in the second part. We start 
with the following observation: If we start at the first index and go right, we see that 
after an index —1 we are always in the second part, while we are always in the first 
part after an index 1. An index zero has no effect. 

First we determine the sign a of the most complicated harmonic sum in M ^^yj^ , nj . 
Given: a harmonic polylogarithm Iljn{n). 

Find: the sign a of the most complicated harmonic sum Sain) in M ^^yg^, • 

We start with cj := — 1 since there is a minus in front of all integrals in the recursions. 
It depends on which part of the recursion we are how the sign is changed: 



Part 1 ((—1)* is present) 




Part 2 ((—1)* is not present) 




(T — >■ cr 
stay 



a 



-a 



a 



a 



a 



-a 



stay goto part 2 goto part 1 stay 



a —7- —(7 
stay 



Analyzing these diagrams we notice that a zero always changes the sign, an other index 
just changes the sign, if the next nonzero index to the left of it is —1. In addition, the 
first nonzero index changes the sign if the factor 1/(1 + x) is present, otherwise if the 
factor 1/(1 — x) is present, the sign remains unchanged. Summarizing, it is enough to 
count the indices less or equal zero. If there are k indices less or equal zero then the 
sign is: 

^ ^ I -(-1)^, if l/(l-x) is present ^^^^ 
1 (-1)^, if l/(l+x) is present. 

Now we determine the index set of the most complicated sum: 



Given: a harmonic polylogarithm Hmin). 

Find: the index set a of the most complicated harmonic sum Sa{n) in M ^ ^^^^ , nj ■ 

Note that the index set of the harmonic polylogarithm and its most complicated har- 
monic sum are closely related if we use the modified notation (Notation 2.1.4) of har- 
monic sums as follows. From the recursions in the Lemmas 3.4.13, 3.4.14 and 3.4.15 
we get that the index set of the harmonic polylogarithm and the index set of the most 
complicated harmonic sum in the Mellin transform of the polylogarithm are almost 
equal. However, some signs may change and there is an extra index. This extra index 
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arises since the recursions result in a harmonic sum of the form Sk{n), /c = ±1 in the 
base case (just 1/(1 — x) or 1/(1 + x) and no harmonic polylogarithm is present). So 
we append the index 1 to the index set. How are the signs going to change? This 
depends again on the part of recursion we have to consider. If we are in the first part 
(of the diagram above) the sign is not changed, if we are in the second part, the sign is 
changed. 

So we can proceed as follows: We go from right to the left in m. If the next nonzero 
index to the left is 1, then the actual sign of the index stays the same, otherwise if 
it is —1, the sign changes. The leftmost nonzero index keeps its sign if the factor is 
1/(1 — x), otherwise if the factor 1/(1 + x) is present, the sign is changed. After this 
we switch the usual notation of harmonic sums. 



Algorithm 2 Computation of the Most Complicated Harmonic Sum in the Mellin 
Transform of a Weighted Harmonic Polylogarithm 



procedure MostComplicated( jq^^) 
a = m 

Append the index 1 to a 
for i = Length(a) to 1 do 
if a[i] then 

j=i-i; 

while j > and a[j] = do 

j=j-i; 

end while 

if j > and a[j] < then 

a[i] = —a[i] 
end if 

if j = and w = 1 then 

a[i] = —a[i] 
end if 
end if 
end for 

k =number of indices < in m 
return w;(— l)*^5a(ra) 
end procedure 



> w = ±1 



Example 3.5.16. Consider Hi -ifl^ifi-i{x) / (1 — x). There are 4 indices less or equal 
1, so the overall sign is —1 by (3.9). Start with {1,-1,0,1,0,-1,1} as the index set 
for the harmonic sum; we added the extra index 1. Looking at the signs of the indices 
we arrive at {1,-1,0,-1,0,-1,-1}. So the most complicated harmonic sum in the 
Mellin transform o/ Hi _i,o,i,o,-i(a;)/(l - x) is -^i _i _2 -2,-i(n). 

Summarizing, we end up at Algorithm 3.5.2.1 to find the most complicated harmonic 
sum and its sign in M ( ^3^' ^) • 
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Remark 3.5.17. Due to our considerations we notice that the mapping which maps a 
weighted polylogarithm onto the most complicated harmonic sum in its Mellin transform 
is infective. 



3.5.2.2 The Reverse Direction 

Wc arc now able to compute the most complicated harmonic sum in the Mellin transform 
of a harmonic polylogarithm. Now we are ready to consider the following problem. 



Given: a harmonic sum Sa{n). 

Find: a harmonic polylogarithm Hm(x) such that Sa{n) is the most complicated har- 



Analyzing Algorithm 3.5.2.1, we notice that the number of negative indices of the most 
complicated harmonic sum in the Mellin transform of H^(a;)/(1 — x) is even no matter 
how m looks like. On the other hand, this number is odd if we consider the Mellin 
transform of H^(a;)/(1 + x), again it does not matter how m looks like. So we can 
decide how the polylogarithm has to be weighted in order to get a special harmonic 
sum, say Sain), to be the most complicated one in the Mellin transform of the weighted 
polylogarithm. Knowing this, we look at Sa{n) in the expanded notation and do the 
reverse of Algorithm 3.5.2.1: First we drop the last index, then we go from left to 
right; the first nonzero index changes its sign if the logarithm is weighted by 1/(1 + x), 
otherwise if it is weighted by 1/(1 — x), the first nonzero index holds its sign. The 
second nonzero index changes its sign if the first nonzero index is negative in the 
current configuration, otherwise it holds its sign. Proceeding in this manner to the last 
nonzero index we will get the index set m, such that Sa{n) is the most complicated 
harmonic sum in the Mellin transform of Hmix) weighted by the already known factor. 
It remains to decide about the sign of the polylogarithm. In order to do that, we can 
use (3.9) where k is the number of indices less or equal zero in m. 

Let us summarize our considerations; compare [33]: Given a multiple harmonic sum 
Sain), wc want to compute a polylogarithm Hm(a;)/(1 ibx) such that Sain) is the most 
complicated harmonic sum in the Mellin transform of H^(a;)/(1 ± x). Let us look at 
the expanded index set of the harmonic sum: 

• If the number of negative indices is even, set a = —1, and there will be a factor 
/ = 1/(1 — x); otherwise set a = 1, there will be the factor / = 1/(1 + x). 




• Drop the last index. Take the remaining indices as the indices of the polyloga- 
rithm. 
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• The first (leftmost) nonzero index clianges its sign if/ = l/(l + x), otherwise if 
/ = 1/(1 — x), the first index holds its sign. 

• Working from the second nonzero index to the right, each nonzero index will 
change its sign if the nonzero index to the left is negative, otherwise it holds its 
sign. 

• Multiply the term by a{—l)^ where k is the number of indices which are less or 
equal zero in the current configuration. 

Remark 3.5.18. Due to our considerations we notice that the mapping, which maps a 
harmonic sum Sain) onto a weighted harmonic polylogarithm Hm(a:)/(1 zizx) such that 
Sa{n) is the most complicated harmonic sum in the Mellin transform of this ilm{x)/{l± 
x), is injective. Hence, using Remark 3.5.17, it follows that the mapping which maps a 
weighted polylogarithm onto the most complicated harmonic sum in its Mellin transform 
is bijective. 

3.5.3 Computation of the Inverse Mellin Transform 

The computation of the inverse Mellin transform of a harmonic sum now is straight- 
forward [33]: 

• Locate the most complicated harmonic sum. 

• Construct the corresponding harmonic polylogarithm. 

• Add it and subtract it. 

• Perform the Mellin transform to the subtracted version. This will cancel the 
original harmonic sum. 

• Repeat the above steps until there are no more harmonic sums. 

• Let c be the remaining constant term and replace it by M.~^(c), or equivalently, 
multiply c hy 6(1 — x) (see the beginning of this section). 

A detailed description is given in Algorithm 3.5.3. 

Example 3.5.19. Consider the sum S'i,i,2,i(ra), which is S'i,i,o,i,i(?T') in the expanded 
index set notation (see Notation 2.1.4). This sum is the most complicated harmonic sum 
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Algorithm 3 Inverse Mellin Transform 

procedure lNvMELLlN(S'a(n), x) > x is argument of output function 

b = a 

i =number of negative indices in 6 
if i is even then 

f = l/{l-x) 

a = -1 
else 

f = l/{l + x) 

(7=1 

6[1] = — &[1] o first index of & changes sign 

end if 

for z = 2 to Length(6) do 

if b[i - 1] < then 
b[{\ = -b[{\ 

end if 
end for 

transform Sh{n) to its expanded notation > See 2.1.4 

drop the last index of b 

k =number of indices < in a 

m = Sa{n) - C7(-1)'=MELLIN(/H„(a;),n) 

while M contains a harmonic sum do 

s =most comphcated harmonic sum in M i> Use 3.5.6 

M = M - InvMellin(s, a;) 
end while 

multiply the constant term in M by 5(1 — x) 
return M 
end procedure 
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in the Mellin transform ofILi^ifi^i{x)/{l—x). The Mellin transform of lli^ifi^i(x) / {1— x) 
yields: 

-4C5 - C25'i,i,i(n) + 5'i,i,2,i(n). 
Therefore after adding and subtracting we arrive at: 

Hi,i,o,i(a;) 



1 -X 



+ 4C5 + C25l,l,lH 



Now we have to consider Si^i^i{n). This sum is the most complicated harmonic sum in 
the Mellin transform of —Hi^i^x) /{I — x). The Mellin transform o/— Hi^i(a;)/(1 — a;) is 
Si,i,iin). We get: 

Hi,i,o,i(x) Hi,i(a;) 

— z \- 4C5 - C2-; • 

1 — x 1 — X 

There is no harmonic sum left. So we multiply the constant term by 5{l — x) and arrive 

at the inverse Mellin transform o/ S'l, 1,2,1 (n) : 

1 — X 1 — X 



With HarmonicSums this can be carried out as foUows: 



HarmonicSums session. 
in[33]- InvMellin[S[l, 1, 2, 1, n], x] 

z2H[l,l,a:] H[l,l,0,l,x] 



Out[33]= 4 z5 Deltalx + 



-1 + x -1+x 



in[34]:= % // ReplaceByKnownFuiictions 

, , , z2H[l,l,a;] H[l, 1, 0, 1, a] 

Out[34]= 4 Deltalx z5H } \ 

— 1 + a; — 1 + a; 

in[35]~ InvMellin[S[2, -1, -1, n], x] 

(-1)"H[-1, 1]^H[0, x] + (-1)"H[-1, 1]H[0, -1,1]- (-1)"H[-1, 1]H[0, -1, x] 

Out[35]= — + 

1+X 

^(H[-l, 1]^H[0, x] - H[-l, 1]H[0, 1, 1] + H[-l, 1]H[0, 1, x] + H[0, 1, -1, 1] - H[0, 1, -1, x]) + 
Deltalx(-H[-1,1]^H[-I,0,l] - H[-l, 1]^H[0, -1, 1] + H[-l, 1]H[1, 1]H[0, 1, 1] 
H[-1,1]2h[1,0,1] + H[-1,11H[-1, 0,-1,1] - H[1,1]H[0, 1,-1,1] - H[-l, 1]H[1, 0, 1, 1] + 
H[1,0, 1,-1,1]) 

in[36]— % // ReplaceByKnownFunctions 

^ , , „ Deltalxln2'* 5 Deltalx ln22 z2 In2z2 (-I)"ln2z2 Deltalxz22 
Out[36i= Deltalxh4half + + + ^-^ + \ ^ + + 

^3 ^ ^^,2 / 1 ^ tlL] HfO x] - (-l)"l'^2H[0,-l,x] ^ ln2H[0,l,a;] ^ H[0,l,-l,x] 
—1 + x \1 — x 1 + x) 1+x 1 — x —1 + x 
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3.6 Differentiation of Multiple Harmonic Sums 



We are now able to calculate the inverse Mellin transform of harmonic sums. It turned 
out that they are usually linear combinations of harmonic polylogarithms weighted by 
the factors 1/(1 it x) and they can be distribution- valued. By computing the Mellin 
transform of a harmonic sum we find in fact an analytic continuation of the sum to 
n G M. For explicitly given analytic continuations see [11, 4, 5, 12, 6]. The differentiation 
of harmonic sums in the physic literature has been considered the first time in [11]. 
Worked out in [6, 8, 7] this allows us to consider differentiation with respect to n, since 
we can differentiate the analytic continuation. Afterwards we may transform back to 
harmonic sums with the Mellin transform. Differentiation turns out to be relatively easy 
if we represent the harmonic sum using its inverse Mellin transform as the following 
example suggests (see [30, p. 81]): 

Example 3.6.1. 



d 

dn 



[ x"H_i(x)dx= / x"log(a;)H_i(x)da; = [ x"H_i,o(x)dx + / x^Hq _i(x)dx. 
Jo Jo Jo ' Jo ' 



Based on this example, if we want to differentiate Sain) with respect to n we can 
proceed as follows: 

• Calculate the inverse Mellin transform of Sa{n) . 

• Set the constants to zero and multiply the remaining terms of the inverse Mellin 
transform by Ho(a;). This is in fact differentiation with respect to n. 

• Calculate the Mellin transform of the multiplied inverse Mellin transform of Sa{n) . 
Example 3.6.2. Let us differentiate S2,i{n). The inverse Mellin transform is: 

5mM^ + M-1(2C3). 

1 — X 



Hence we have: 

S2,l{n) = /'^^n Ho,l(^)-C2 ^^^ / ^nM-l(2C3)dx 



r --Ho,(x)-Ho,(l)^° ^ ^^^^ ^ ."M-(2C3)cZ: 
Jo 1 — X Jq 1 — X Jq 



M|1^.„)-m(^,„,.2C3. 
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Differentiating the right hand side with respect to n yields: 



Jo 1 - a; Jo l-x 

= [\x-^^^+X-^^^-C2X-^dx 

Jq 1 — X 1 — X 1 — X 



^ /"^^n HoAil^) -Ho,o,i(l) I ^n Ho,i,o(x) -Ho,i,o(l) ^^^, Ho(x) -Ho(l) ^^ 
Jo 1 — X 1 — X 1 — X 

1 ^riTT f„\ TT /I \ nl 



2 , x"Ho,o,i(x)-Ho,o,i(l) ^^ ^ p x"Ho,i,o(x)-Ho,i,o(l) ^ 



,rx-Ho,o,(x)-Ho,o,(l)^^^ /■ 
Jo 1 - a; Jo 

-cJ'-^""°W-"°Mdx 



^ + C2S2(n)-S2,2(n)-2S3,,(n). 



Algorithm 4 Differentiation of harmonic sums 

procedure DlFFERENTIATEHARMONIcSUM(5a(n)) 

h = lNVMELLIN(S'a(n), x) 

replace 5{l — x) in /i by 
h = Ho(x) * h 

expand all products in h > Use (3.3) 

nsum =number of summands in h 
for i = 1 to nsum do 

replace iih summand of h by MELLiN(ith summand,ra) 
return h 
end for 
end procedure 



The package HarmonicSums provides a tool to differentiate harmonic sums: 



HarmonicSums session. 




in[37]:= DifFerentiateSSum[S [2, 1, 


n] , n] // ReplaceByKnownFunctions 


Out[37i= ^+z2S[2,n]-S[2,2,n] 


- 2S[3,l,n] 


in[38l~ DifferentiateSSum[S[5, n] 


n] // ReplaceByKnownFunctions 


8 z2^ 

Out[38]= — 5S[6,n] 




in[39]:= DifferentiateSSum[S [2, -1 


2, n], n] // ReplaceByKnownFunctions // Expand 
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,r „ In2*z2 , o „2 121 z2^ 11 z3^ , ln2*S[2,n] 

Out[39]= 4s6 + 41i4halfz2 H In2^z2^ — 41i4halfS 2,n + 

^ 6 120 4 6 

ln2^z2S[2,n] + "^^ ^^^q ""^ + 2z3S[2,-l,n] - S[2,-2,2,n] - 2S[2,-l,3,n] - 2S[3,-l,2,n] 



3.6.1 Application of Differentiation 



As shown in [6, 8, 7] new relations between multiple harmonic sums arise if we introduce 
differentiation on harmonic sums. We can now use these relations to reduce the number 
of basic sums we computed in the first strategy of Section 2.7 by allowing relations due 
to differentiation. Unfortunately we cannot adopt the second strategy since we have to 
know the actual values of the indices of a harmonic sum in order to differentiate it. Let 
us look at an example. 

Example 3.6.3. Using only algebraic relations we have to take Si{n) and S2{n) into 
the basis, but using in addition the differential operator ^ we get the following relation: 

So we can express 5*2 (n) with the help of Si{n) and we can eliminate 5*2 (n) from the 
basis. In general we can express all sums of depth one by ^S±i{n) with i > 0; see [6j. 



As indicated in the previous example for S2{n) we proceed as follows. First we fix 
an algebraic basis in which we can represent our expression (see the first strategy in 
Section 2.7). For each multiple harmonic sum Sa{n) in the algebraic basis we try to 
find another multiple harmonic sum Sb{n) such that Sb{n) is less complicated than 
Sa{n) and such that Sa{n) is a most complicated harmonic sum in ^Sf,{n). We can 
do this by computing the inverse Mellin transform of Sa{n) = 501,02, .■•,afc('^)- There 
we look for the most complicated harmonic polylogarithm. We call a harmonic poly- 
logarithm most complicated if it is the polylogarithm with the largest weight or in the 
case of identical weights it has the largest number of nonzero indices [33]. If at least 
one of the indices of the most complicated harmonic polylogarithm Hm(-T) is zero (so if 
Hm(a^) = ^0,7712(2^))) then we can find a harmonic sum Sb{n) with the desired prop- 
erty, since wc differentiate by multiplying Ho(x). The harmonic sum Sb{n) is just the 
most complicated one in the Mellin transform of }lmi,m2 i^-) weighted by the same factor 
as Umix) in the inverse Mellin transform of Sain). Hence, due to the consideration in 
Section 3.5 it is clear that the found harmonic sum satisfies the desired property. In 
fact, if there is an index a; of Sa{n) which is not equal to ±1, we will find Sf,{n) with 
6= (ai,...,,a;_i,sign(a/)(|a/| - 1), a^+i, . . . , afe). 
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Weight 


Number of 


Sums 


a-basic sums 


d-basic sums 


1 


2 


2 


2 


2 


6 


3 


1 


3 


18 


8 


5 


4 


54 


18 


10 


5 


162 


48 


30 


6 


486 


116 


68 



Table 3.1: Number of harmonic sums, and the respective numbers of basic sums by 
which all sums can be expressed using the algebraic (a-basic) or algebraic and 
differential relations (d-basic) in dependence on their weight; see [5, 6, 8, 7]. 

Example 3.6.4. We consider the sum S4^2{'n). We can use both -S'3,2(?^) or S'4,i(n) to 
get relations for 84^2 (n) '■ 

J^53,2(n) = 2C| + 2C353(n)-253,i(n)-354,2(r'^ 



= ^-Cl + C254(n)-54,2(n)-4S5,i(n) 



We use all these new relations to reduce the number of sums in the basis as far as pos- 
sible. Summarizing, Si,{n) is less complicated than Sa{n) (which wc want to eliminate) 
and all additionally introduced sums arc not more complicated than Sain). In this way 
wc can successively eliminate the original given algebraic basis to a basis whose sums 
arc algebraically independent and where as much sums as possible can be eliminated 
by the differential operator. In Table 3.1 we can see how the number of basic sums 
reduces in comparison to the use of algebraic relations only; compare [8, 7]. For the 
package HarmonicSums tables up to weight 6 are available. 



HarmonicSums session. 

in[4o]:= ReduceToBasis[S[l, 2, -2, n], UseDifferentiation True] 

7z2z3 5z5 2z2Ditf[.S[-L ;;]w;,.2] DifitS[-i. w;, i] Ditt'[S[-i. ;;],;(, 2; Ditt]S[L /;]. n, 1] 



Out[40] = 



12 24 
z2 Diff[S[l, n], n, 2] 

12 



3 72 

Diff[S[-l, n],n, 1] Diff[S[l, n], n, 2] 

6 



3 

Difr[S[-2. l,nl. n,2l 



piff[S[2,-^^ 
S[-2,2,l,n] 



Here DiS[f{x),x,i] := 
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Chapter 4 



Half-Integer Relations 



In this chapter we allow besides the upper index n, the index | (or equivalently the index 
2n) in the multiple harmonic sums. As it turns out, this leads to additional relations 
by using the operator Half[/(a;), x] = /(§)• Such relations have been determined in 
[6, 8, 7, 11] by exploiting properties of the Mcllin transform. Here we follow a different 
approach. We start with the following lemma, it has already been used in [11]. 

Lemma 4.1.1. Let n,m,a E N. Then the following relation holds: 

5„(2n) - S_„(2n) = ^Sa(n). (4.1) 



Proof. For n G N, 

2n 



5_„(2n) - g - g - (2^^i)a) - ^-^"H - E (2i-l)a 

2"^'^«H ~ ((^ + (2i - l)a ~ (2i)^) 



^S'„(n) - S'a(2ra) + ^S'a(n) 



-^5„(n)-5„(2n). 



□ 



Theorem 4.1.2. (see /-^-Z/y) i^ei n, m G N and G N /or z G N. T/ien w;e have the 
following relation: 

'S'±am,±a^_i,...,±ai (2n) = ^-^^^ ai-m ^'^m.,a.m-i,—,a.i l'^) ('^•2) 

where we sum on the left hand side over the 2"* possible combinations concerning it. 
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Proof. Wc proceed by induction on the depth m of the multiple harmonic sums. In 
Lemma 4.1.1 we have already proven the case m = 1. Let's assume (4.2) holds for m: 

" / 1 



n / / -,x2i 



1+ 

n 



_ ]^)a„+i X] 'S'±a„,±a^_i,...,±ai(2n)J 



= 2 5^^^^^^S±a^,±a^_i,...,±ai(2n) 
" 1 1 

^ 2EII1 ai-(m+l)'^'^'"+i'«""-'«i*^"'^- 

Hence (4.2) holds for m + 1. □ 
Example 4.1.3. For ra G N, 

ls2,3{n) = 52,3(2n) + S2,-3(2n) + 5-2,3(2n) + 5_2,-3(2n). 



Similar to Section 3.6.1 we use all these new relations to reduce the number of sums in 
the basis as far as possible. In the Tables 4.1 and 4.2 we can sec how the number of basic 
sums reduces if wc use algebraic relations and half-intcgcr relations in comparison to the 
use of algebraic relations only. If wc use all three types of relations, namely algebraic, 
differential and half-integer relations, the number of basic sums reduces further as we 
can see in the Tables 4.1 and 4.2. 

Example 4.1.4. From Example 3.6.3 we know that all harmonic sums of depth one 
can he expressed by the two sums Si{n) and S-i{n) by using differentiation. From 
Theorem we get: 

5_i(2n) = 5i(2n) -5i(n). 

Hence we can express S-i{n) by Si{n) and S\{n/2), and we can represent all harmonic 
sums of depth one by Si (n) using differentiation and half-integer relations. 



For the package HarmonicSums tables using algebraic and half-integer relations up to 
weight 6 are available. We remark that all sums (d-, h-, dh-basic sums) in the Tables 
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Number of 


Weight 


Sums 


a-basic sums 


d-basic sums 


h-basic 


sums 


dh-basic sums 


1 


2 


2 


2 




1 


1 


2 


6 


3 


1 




2 


1 


3 


18 


8 


5 




6 


4 


4 


54 


18 


10 




15 


9 


5 


162 


48 


30 




42 


27 


6 


486 


116 


68 




107 


65 



Table 4.1: Number of harmonic sums, and the respective numbers of basic sums by 

which all sums can be expressed using algebraic (a-basic sums), differential 
(d-basic sums) and/or half-integer relations (h-basic/ dh-sums); see [5, 6, 8, 
7]. 





Number of 


Weight 


Sums -1 { — 1} 


a-basic sums 


d-basic 


sums 


dh-basic sums 


1 


1 


1 




1 


1 


2 


3 


2 




1 





3 


7 


4 




2 


2 


4 


17 


7 




3 


3 


5 


41 


16 




10 


9 


6 


99 


30 




17 


17 



Table 4.2: Number of sums, which do not contain the index {—1} and the respective 
numbers of basic sums by which all sums can be expressed; see [5, 6, 8, 7]. 



4.1 and 4.2 are algebraically independent by construction. In [8] the representations of 
the dh-basic sums (up to weight 5) have been established with a different approach. 

HarmonicSums session. 

in[4ii:= ReduceToBasis[S[l, 2, -2, n], UseHalflnteger True] 

o.[«H S[l,n] -S[2,n]) S[2,n] + S[l,n] - S[4,n]) - S[-4,l,n] - 

gj^^^j j- Half[S[4, n],n] ^ j^ Half[S[2, n], n] _ g^^^ S[2, n] - S[4, n] - S[2, -2, n]) 
^ Half[S[2, n],n] _ g^^, s[2, 1, n] + S[3, -2, n] + S[-2, 2, 1, n] 



Moreover tables using all three types of relations up to weight 6 are available. 



HarmonicSums session. 

in[42]:= ReduccToBasis [S [1 , 2, -2, n], UseDifferentiation True, UseHalflnteger True] 

7z2z3 5z5 2z2I)iff[lla.lf[S;i. ;ij - S;i, /(. 2] Diff[llalf[S[i, ;i] - S^i, /i, 4) 
Out[42]= — — 24 3 72 ~'~ 

Diff[Half[S[l, n],n] - S[l, n], n, 2] Diff[S[l, n],n, 1] _ z2 Diff[S[l, n], n, 2] _ 
3 12 
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Diff[Half[S[l,n],n] -S[l,n],n,l]Diff[S[l,n],n,2] Diff[S[-2, 1, n], n, 2] Diff[S[2, -2, n], n, 1] 
6 6 3 

+ S[l,n]S[2,-2,n] + ^^^^^'^'"^ + Diff[Half[S[l, n], a] - S[l, n], n, 1] S[2, 1, n] + 

S[-2,2,l,n] 



Further investigations in this direction in connection with Chapter 6 are in 
progress [1]. 
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Chapter 5 

Summation of Multiple Harmonic 
Sums 

With the package Sigma [36] we can simphfy nested sums such that the nested depth 
is optimal and the degree of the denominator is minimal. This is possible due to a 
refined summation theory [39, 37] of HS-ficlds [22]. In the following wc want to find 
sum representations of such simplified nested sums in terms of harmonic sums as much 
as it is possible. Inspired by [35] for harmonic numbers and [28, 41] we consider sums 
of the form 

n 

^sr(z) (5.1) 

i=l 

where s G S{i) or s G 1)*] and r{i) is a rational function in i. We can use the 

quasi-shuffle algebra property of the multiple harmonic sums to split such sums into 
sums of the form 

n 

J2r{i)Sa(i) (5.2) 

i=l 

or 

n 

^(-l)V(i)5„(i) (5.3) 

i=l 

where r(i) is a rational function in i. 

In this chapter we will present new formulas to rewrite sums of the form (5.2) and (5.3) 
and we will show how we can use these formulas in combination with Sigma. In the 
following we consider the problem: 
Given: a sum a of the form (5.1). 

Find: as much as possible a representation of a in terms of harmonic sums. 
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5.1 Polynomials in the Summand 



First we consider the case that r{i) is a polynomial in Z[i], i.e., r{i) = Pm'i'^+- ■ ■+Pii+Po 
with pk E 7.. Here we generalize the harmonic sum case from depth 1 [35] to arbitrary 
depth. If wc are able to work out the sum for any power of i times a multiple harmonic 
sum (i.e., Y17=i ^ aii) or ^iLi(^l)*^™5'a(^)5 € N), wc can work out the sum for 
each polynomial r{i). Let us start with m = where the factor (—1)' is present. 

Theorem 5.1.1. LetneN and au G Z/{0}. Then 

J2 i-^ySa„a2,...{i) = \ (5-ai,a„...(n) + {-Xf Sa„a2^...{n)) . (5.4) 



Proof. For n & N 

n 



i=l 



i=i j=i 
sign(ai'' 



E^s«...-.o)i:(-i)' 



n 



sign^ai 



1 



J 



i-iy + i-iy 



(5-a„a„...(n) + (-l)"5„„„„...(n)). 



□ 



When m = and the factor ( — 1)* is not present wc have to distinguish between several 
cases. In order to find this identities C. Schneider's Sigma package played a decisive 
role. 

Theorem 5.1.2. Let k,n,l G N and aj G Z/{0}. Then Y17=i '^"■i,a2,..;ai{i) can be 
simplified as follows. 



\ai\ > 2 : 



a 



0-1,0,2, 



„a,{i) = {n + l)S, 



01,02, 



■,oi 



(n) - 5'sign(ai)(|ai|-l),a2,...,a(H; 



ai = — 1 : 

n ^ 

'5-1,02,03,. ..,0; (^) = (^+l)'S'-l,a2,a3,...,ai (^^) — 2 {^-a2,ai,...,ai{'n) + ( — I)"5a2,a3,...,ai (^)) ] 

1=1 
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aj = 1, I < j < I ■ 

n 

1] 5i, 1, . . , , i(i) = (n + 1) (^1, 1, . ._^(n) - 5i, 1, . ._^(n) + • • • + (-l)'-^5i(n) + (-1)') 



I times l — l times 

-(-ly+i; 



aj = 1, l<j<l, ai+i = -1 : 

n 

E^M^,-i» = (n + l)(5i^i_^,_i(n)-5i^._^,_i(n) + --- + (-iy5_^ 

^ i times I times l — l times 

(-1)' 

+^(i-(-ir); 

aj = 1, l<j< l,ai+i = -1 : 

n 

I times l — l times 



a-j = 1, 1 < j < ^, |a/+i| > 2 : 

n 

Xl'S'l,l,-.-,l,ai+i,«!+2,-W = (" + l)('S'l,l,...,l„a,+i,ai+2,...H " -^l, l, . . . , l,a,+i,a(+2,... (") 



I times l — l times 

+ (-l)'+^Ssign(a(+i)(l»m|-l),«i+2,-(^)- 



Proof. We give a proof of the first and tlie third identity all the other identities follow 
similarly. First let |ai| > 2: 



i=l j=l •> 



7>i| 

i=l i=l j=l 



7>i| 



^ sign(ai)-' / -N/ , 1 •^ 

j=l J 

= (n + l)Saj_,a2,...{n) - 5'sign(ai)(|ai|_i),a2,...H- 
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Consider the third identity, i.e., aj = 1, 1 < j < I'- We proceed by induction on the 
depth I. For 1 = 1, the right hand side is (n + l)Si{n) — n and the left hand side gives: 



E^i« = EE7 = E7Ei = E>+i-^-) = (-+i)Mn)-Ei 

i=l i=l 3=1 •' j=l •' i=j j=l •' j=l 

= {n + l)Si{n) -n. 
Now assume that the third identity holds for depth < / — 1. Then 

n n i ^ " 1 " 

l-l times l-l times 

n ^ n 

= E /m^(j) (n + 1 - i) = (n + l)Si^i_i (n) - ^ 5i_i, . . . , i (i) 

l-l times i times I- 

= (ri + l)Si,i,...,iH-(n + l)(5i,i,...,i(n)-5i,i,...,i(n) + 



i — 2 times 



+(-iy-25i(n) + (-1)^-1) -(-1)' 

(n + l)(5i,i,...,i(n)-5i,i,...,i(n) + --- + (-iy-iSi(n) + (-iy) 



/ times l — l times 



+(-1)' 

Therefore the identity holds for I and we finished the proof of the third identity. □ 



In general, for m > 1 we can use the following theorem. 
Theorem 5.1.3. Let k, neN and ai G Z/{0}. Then 

1=1 i=l j=l i=l 

n n n ^ , j-1 

Ei-^y^'Sa,,a,,M = Sa„a,,..Xn) ^ (-l)'^^ - ^ ^"^'^Z-^'^ ^ (-1)^^^ 
i=l 1=1 j=l •' i=l 



Proof. We just give a proof for the first equality. The second follows analogously. 

V^.fcr, /-N _ ,k \^ sign(ai)J5'aa,a3,...(0 
1=1 1=1 j=l ■' 

sign(aiV.ga,,„3,...(i) A 

Z-^ j\ai\ 

j=l i=j 
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sign(aip5a,,a3,...(0 / ^ -fc _ ^ -fc 



J 



i=l j=l i=l 



□ 



Since the sums J27=i ("1)*^*^ s^nd can be expressed as a polynomial in n, to- 

gether with a factor (—1)", we can apply Theorem 5.1.3 recursively and we end up 
at sums over harmonic sums that can be handled in Theorem 5.1.1 or Theorem 5.1.2. 
Summarizing, we can always work out the sums Y17=i ^(^)'S'a(O and Y17=i (~l)*''(0'S'o(0 
where r(i) is a polynomial; the result will be a combination of multiple harmonic sums 
in the upper summation index n, rational functions in n and the factor (—1)". 

Example 5.1.4. For n G N 

7n + n^ (3 + 2n-n^) Sijn) _ Si,i{n) 
^ 6 6 

n (1 + n) (l + 2n) S2,i{n) 
6 



E^'^2,i(o = 



i=l 



Our method is built in inHarmonicSums with th function call TransformToSSums. 



HarmonicSums session. 

n 

in[43l- TransformToSSums[^i'' S[2, 1, i]] 

Out[43]= =^(-44n + 81n - 46n^ + 9n^) + (-5n + 2n^ + 4n^ - 3n*)S[l, n] + -^S[l, 1, n] + + 

/ ZU uU oU oU 

n)(l + 2n)(-l + 3n + 3n^)S[2, 1, n] 



5.2 Special Rational Functions in the Summand 



Next we consider the case r{i) = + cf , where c G N is a fixed number. First we 
look at sums of the form Sa{n + c), compare [28]. 

Lemma 5.2.1. Let c,n,k e N, aj G Z/{0} for j e {1,2, . . . ,k} . Then forn>0 



Summation of Multiple Harmonic Sums 



83 



and n > c, 



Sa\,a2,...,au{'^ 



n 



Proof. We only proof the first identity: 

n+c . , 

c _ y^ sign(ai)^ 

'Jai,a2,...,afcl,''' t-; — ^ |^^| <:>a2,...,afc U; 



Q C^^^ sign(ai)^ ^ 

'-'ai,a2,...,afc W -r ^ |^^| iJa2,.--,afc 

j=n+l ■' 

( \ sign(ai)-'''''^ 



□ 



Since the multiple harmonic sum on the right sides of (5.5) and (5.6) have reduced 
depth, we can relate a sum of the form S'a(nib c) to sums of the form Sa{n) by recursive 
application of (5.5) or (5.6). Namely, we can use this lemma to work out sums of the 
form 

y- Sa{i) i Sgji) 

^ (i + c)"» ' (i + c)™' 

We shift the summation index: 

Sa{i - C) i Sa{i-c) 
j=c+l i=c+l 

and afterwards we use (5.6) to relate the sum Sa{i — c) to Sa{i)- Similarly we can work 
out sums with ' — ' in the denominator i.e., 

y or y (-i)^-^« 



i=l+c ^ ' «=l+c ^ ' 

Example 5.2.2. For G inN, 

i=\ 1=2 i=2 

n+1 



+ 52,3(0 

i 



= E 



(-1)^^3(0 (-1)^^2,3(0 



Z3 
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1=2 1=2 
= 5_3,3(n + 1) + 1 - 5_3,2,3(« + 1) - 1 

= S-aM + ^-^p^J^ S_i,2,3(n). 



5.3 General Rational Functions in the Summand 



If we want to work out sums of the form (5.2) or (5.3) for a general rational function 
r(n) = ll^, where p(n) and q{n) are polynomials, we can combine these strategies. 
We start as follows: 



If the degree of p is greater than the degree of q we compute polynomials f(n) and p{n) 
such that r(n) = f{n) + and the degree of p{n) is smaller than the degree of q. We 
split the sum into two parts, i.e., into 

n n 71 _/ .N 

^r(i)5„(i) = Y.f{i)Sa{i) + jJ.Sa{i) 

i=l i=l i=l ^^^> 

or 

II, n n —/■\ 

^ (-l)V(i)5„(i) = ^ (-l)V(z)5„(i) + ^ (-l)^^5„(i). 
j=i i=i i=i 

The first sum can be done using Theorems 5.1.1, 5.1.2 and 5.1.3. For the second sum 
we proceed as follows: 



1. Factorize the denominator q{n) over Q. 



2. Let A be the product of all factors of the form (i + c)"^ with m G N and c G Z, 
and let B be the product of all remaining factors such that A{n)B{n) = q{n). 



3. 



For example with the extended Euclidean algorithm (see Remark 5.3.2) we com- 
pute polynomials s(n) and t{n) such that s{n)A{n) + t{n)B{n) = p{n). This is 
always possible since A and B are relatively prime. Hence we get 
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4. Split the sum into two sums, each sum over one fraction, i.e., 



i=l i=l 



or 



i=i ^^^> i=i ^ i=i ^ 

5. The sum with the denominator B remains untouched (for details see Remark 
5.3.1). We can now do a complete partial fraction decomposition to the first 
summand and split the sum such that we sum over each fraction separately. 



6. Each of these new sums can be expressed in terms of harmonic sums following 
Subsection 5.1. 

7. We end up in a combination of rational functions in n, harmonic sums with upper 
index n, the factor (—1)" and perhaps a sum over the fraction with denominator 
B. 

Remark 5.3.1. In our implementation the sum %Sa{i) is passed further to 

Schneider's Sigma package. The underlying difference field and difference ring algo- 
rithms [39, 38, 37] can simplify those sums further to sum expressions where the de- 
nominator has minimal degree. The result of Sigma is again passed to the package 
HarmonicSums. If possible, it finds a closed form in terms of harmonic sums. 

Remark 5.3.2. In our implementation we do not use the extended Euclidean algorithm 
since it turned out that it is faster to construct a partial fraction decomposition with 
unknown coefficients. After clearing the denominators we have to solve the linear system 
of equations which we get by coefficient comparison. 

Example 5.3.3. Let us give several examples using the package HarmonicSums; the 

first sum is without a 'had' part: 



HarmonicSums session. 

in[44]:= TransformToSSums [ V ^^t^'^''] ] 

3n 3S[2,n] (-11 - n + 3n^) S[3, n] 7S[4,nl (-19n - 20n^ - 5n^) S[2,3,n] 

Out[44]= . r 1 . , , , 7 1 -. 1 



4(l + n) 4 4(l + n)(2 + n) 4 2 (1 + n) (2 + n) (3 + n) 

S[3,3,n] +S[l,2,3,n] 



in[45].= TransformToSSums[V ' ^' ^ 1 
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— + ^ W^) + 

71=1 

3 n S[2, 3, n] + (Eri=i s+br) ^l^, 3, n] + 9 n S[2, 3, n]) + 3 S[3, 3, n] + 3 n S[3, 3, n] 

3 (1 + n) 

3S[l,2,3,n] +3nS[l,2,3,n] 
3 (1 + n) 

Of course we can also work out more difficult nested sums: 

j 1^1 = 1 2^ 

in[46]:= TransformToSSums[y^ 1 

45n + 67n'^ + 66n-^ + 48n'^ + 19n^ +3nP 3 (l + 6n + 4/1^ + n^) S[l,n] , 

Out[46]= F + — ^ 1 — - 2S 2, n + 

4(l + n)" (2 + n) 2(l + n)* ^ ^ 

(-l-6n-4n^-n^) S[l,l,n] (-7 + 3n + 11 + 5n^) S[2,l,n] 9S[3,l,n] 

+ ^ 2(l + nf ^ 2 + 

2 (3 + 3n + n ) S[l,l,l,n] _ 3S[2,l,l,n] - 2S[2,2,l,n] - 3S[3,l,l,n] - 2S[3,2,l,n] + 
(1 +n) 

2S[2,l,l,l,n] +2S[3,l,l,l,nl 



We can also work out S-sums: 

in(47l- TransformToSSumsfy^ - ^ 1 

-3 (2-2^+"3"+4n-6"n + 2n') 3^+" (7+10n + 4n^) S[l,{2},n] _ 5S[l,{6},n] 
{1 + nf ^ (2 + 3n + n2)2 2 

^^^'^'""^ - S[3, {6}, n] + S[2, 1, {3, 2}, n] 



Out[47] = 



Of course not all sums can be rewritten in terms of harmonic sums: 

n 1 

ini48l:= TransformToSSums ] 



Out[48]= V ( -TTTT^ 



+ 1) i!i 
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Chapter 6 



An Example from Particle Physics 



Single scale quantities as anomalous dimensions and hard seattering cross section in 
renornializable quantum field theories are found to obey difference equations of finite 
order in Mellin space. In [9, 10] a formalism is established to construct general solutions 
for these quantities from a sufficiently large number of fixed moments. From the given 
moments a recurrence relation is established and solved in terms of the Mellin parameter 
N, finite harmonic sums and their generalizations at finite values of N and in the limit 
N — > GO. In these computation our package HarmonicSimis contributed as follows. 
E.g., for the C^C^TVp-term of the 3-loop non-singlet splitting function P^s2i t^]' 
Sigma solved a recurrence of order 7 and found the following closed form evaluation [9, 
Exp. 5]. 

^ •= ^NS,2,CaCfNp = 

2 (l086Af'^ + 3258A^^ + 2129iV5 - 288iV^ - 67N^ - 206N^ - 156iV + 144) 

277V4(Ar + 1)3 

32 {8N^ + 33iV3 + 53Af2 + 25N + 3) ^ 16 A 1 32 ^ (-1)* 

m(N + 1)4 ^~ ^ +yZ^^ + yz^~^ 

i=l i=l 

_ 16 (l0iV2 + loiV + 3) A (-1)^ 1336 1 _ 64 {SN^ + 8N + 3) ^ {-If 



9N{N + 1) ^ i3 '27-^^2 



i=l \ > J 



16 (4Ar6 + 88N^ + 3UN^ + AUN^ + 201 A^^ + iqj^ _ ^2) ^ (_i)i 
^ 9iV2(iV + l)2(Ar + 2)2 

( 8(l4Ar2 + 14jV + 3) l6J^l\J^l 64 ^ F .g^ ^W^ 

3MiV + l) 3^J^i3+3Z^ i +-^2^ (^ + 2)2 

\ ^ ' 1=1 / 1=1 1=1 1=1 ^ ' 



32 {22N^ + 22Ar - 3) ^ E]=i / 32 (2jV2 + 4jV + l) 

m{N + l) ^ i + 2 ^[t^ijx 3(Ar + l)3 ^ 



An Example from Particle Physics 



88 



128 ^(-1)* 32(2iV3 + 2Ar2_3iV_2) ^ (-1)^ 64 A E5=i ^ 



^ i 3N{N + 1){N + 2) ^ i'^ 3 ^ i + 2 / 

256f (-i)'£-..l , mf (E-..y)(EL}) 

3 ^ i 3 ^ i + 2 ■ ^ ■ ^ 

1=1 1=1 



Given such sum representations, the package Harmonic Sums provides several sophisti- 
cated functions to simpUfy such expressions. First wc rewrite (6.1) into an expression 
in terms of harmonic sums (harmonic sums arc defined in Chapter 2) by using the 
function call S :=TransformToSSums[P] (for details on this function see Chapter 5): 

^ = ^NS,2,CaCfNp = 

-2 (144 - 12 n - 362 - 417 - 1603 - 96 (-1)" - 1855 - 384 (-1)" n^) 

27n4(l + n)'' 
-2 (347 + 1080 + 270 n^) 64S'_4(n) 
27n4(l + n)^ ^ 3 

16 S-3{n) (-18 - 13 n + 54 + 34 + 36 n Si (n) + 54 Si (n) + 18 Si (n)) 

9n (1 +n) (2 + n) 

16 5'_2(n) (-6 + 17n + 42n2 + 16n3 - 12nSi{n) - 42n^ Si{n) -6n^Si{n)+4) 



9n2(l + n)^ (2 + n) 



16 5-2 (n) {A8n^Si{n) + 2An^Si{n)) 133652(n) 8 (3 + 14n + 14n=^) 53(n) 
^ 9n2(l + n)2 (2 + n) ^ 27 3n (1 + n) 

16S4(n) 1285_3,i(n) 128 (-1 + n + n^) S_2,i(^) 1285i_3(n) 



3 3 3 (1 + n) (2 + n) 3 

32 (-6 + 5n + 42n2 + 22n=^) Si_2(n) 4Si(n) (18 + 36n + 36 n^ + 281 n^) 

9n (1 + n) (2 + n) 27n3(l + n)^ 

45i(n) (72(-l)"n3 + 627 + 627 + 209 n^ + 36n3 53(n) + 108n^53(n)) 

27 n3 (1 +n)-'^ 

45i(n) (108 n^ 53(n) + 36 n^ 53(n) + 144 5i _2(n) + 432 5i _2(n)) 

27n3(l + n)^ 

4 5i(n) (432n5 5i,_2(n) + 144 n^ 5i,-2(n)) ^ 64 5i,3(n) _ 32 52,-2 (n) 



27n3(l + n)'^ 

^ 128 5i,-2,i(n) ^ 128 5i,i,-2(n) _ ^^ 2) 



In (6.2) 15 different harmonic sums appear: 



5_4(n), 5_3(n), 5_2(n), 5i(n), 52(n), 53(n), 54(n), 5_3,i(n), 5_2,i(n), 5i,_3(n), 
5i _2(n), 5i,3(n), ^2 -2(n), 5i _2,i(n), 5i,i _2(n). 
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Wc can use the algebraic relations between harmonic sums derived in Chapter 2 to 
reduce the number of different harmonic sums appearing in (6.2). Namely, activating 
the call ReduceToBasis[5'] we get: 



■ NS,2,CaCfNf ~ 

2 

J (-144 + 12 n + 362 + 417 + 1603 n*^ + 96 (-1)" n*^ + 1855 

27n4(l + n)^^ 

+384 (-1)" - 347 - 1080 - 270 + (288 n'^ + 1152 + 1728 + 1152 rJ 

+288 n^) S_4(n) + S-3{n) (-72 - 456 - 1176 - 1512 - 960 - 240 

+ (144 n'^ + 576 + 864 + 576 + 144 n^) Si (n)) + (668 + 2672 

+4008 + 2672 rJ + 668 n^) S2{n) + 5_2(n) (-72 + 96 + 792 + 1008 

+384 + (-480 n'' - 1920 - 2880 - 1920 - 480 n^) Si{n) + (288 n'' + 1152 

+1728n^ + 1152n^ + 288n^) S2{n)) + (-108 - 828 - 2340 - 3132 - 2016 

-504 n'^) Ss{n) + (360 + 1440 + 2160 + 1440 n'^ + 360 n^) S^n) 

+ (- 144 + 48 + 1488 + 2736 n^ + 1920 + 480 n^) 5_2,i(n) + S'i(n) (-36 

-lOSn^ - 144 - 634 - 144(-l)"n^ - 1816 - 144(-l)"n^ - 2508 

-1672n^ -418n® + (21671"^ + 864 + 1296 + 864 + 216 n®) 83(71) + (288n^ 

+1152 + 1728 + 1152 + 288 n^) S'-2,i(n)) + (-144n^ - 576n^ - 864n6 

-576 ra^ - 144 n^) ^2 -2(n) + (-288 n"^ - 1152 - 1728 ra^ - 1152 - 288 n®) 53,1(71) 

+ (-576 - 2304 - 3456 - 2304 - 576 n^) S-2,1,1 (n)) (6.3) 



In (6.3) only the following 11 different harmonic sums appear: 

S-4in), S-3{n), 8-2(71), Si(n), 82(71), 83(71), 84(71), 8-2,1(71), 82,-2(n), 83^1(71), S_2,i,i(n). 

We want to emphasize that these sums are algebraic independent. Exactly these repre- 
sentations were used in [9] to establish compactified representations originally computed 
by [29] and [42]. Using relations originating from differentiation (for details see Chap- 
ter 3) we can go further. We use the function call ReduceToBasis[5', UseDifferentiation 
True]: 



^NS,2,CaCfNp 



J (720 - 60 n - 1810 - 2085 - 8015 71'^ - 480 (-1)" - 9275 n 

135n4(l + ^^)^ ^ ^ ^ 

-1920 (-1)"* + 1735 + 5400 + 1350 - 180 C2 + 240 (2 - 1360 n"^ (2 
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-10840 C2 - 19080 (2 - 13360 C2 - 3340 C2 + 1224 n'^ (2^ + 4896 (2^ 

+7344 C2^ + 4896 n'^ (2^ + 1224 (2^ + 270 Cs + 2430 n'^ Cs + 7290 Ca 

+9990 ra^ Cs + 6480 ra^ Ca + 1620 Cs + (I8O + 1140 + 2940 + 3780 + 2400 ra^ 

+ (3340 n"^ + 13360 + 20040 + 13360 + 3340 - 720 C2 - 2880 C2 - 4320 

-2880 n^C2 - 720 n^s) B[Si{n),n,l] + (-360 n'' - 1440 - 2160 

-1440n^-360n^) D[5i(n),n, 1]^ + D[S_i(n),n, 1] (-360 + 480 + 3960 n'' 

+5040 + 1920 + 1440 C2 + 5760 (2 + 8640 C2 + 5760 C2 + 1440 C2 

+ (-1440 n'' - 5760 - 8640 - 5760 - 1440 n^) D[5i (n),n, 1]) + (270 

+2070 + 5850 n^ + 7830 n^ + 5040 n^ + 1260 n^) D[5i(n),n,2] + (360 + 1440 

+2160 n^ + 1440 n^ + 360 n^) D[5i(n),n,3] + (-720 - 2880 - 4320 - 2880 

-720 n^) D[S'2,i(n),n, 1] + (720 - 240 n"^ - 7440 - 13680 - 9600 - 2400 n^) 

5-2,1 (n) + Si{n) (180 n + 540 + 720 n''^ + 3170 + 720 (-l)"n^ + 9080 

+720 (-1)" + 12540 + 8360 n'^ + 2090 - 1200 C2 - 4800 (2 - 7200 C2 

-4800 C2 - 1200 (2 - 540 ri^ Cs - 2160 Cs - 3240 Cs - 2160 Cs 

-540n^Cs + (-2400 n'' - 9600 - 14400 - 9600 71"^ - 2400 n^) D[5_i (n),n, 1] 

+ (-360n^ - 1440n^ - 2160n^ - 1440n^ - 360n^) D[5_i(n), n, 2] + (-540n^ - 2160ra^ 

-3240 - 2160 -540 n^) B[Si{n),n,2] + (-1440 - 5760 - 8640 - 5760 

-1440 n^) 5_2,i(n)) + (720 + 2880 + 4320 + 2880 + 720 n^) S2 -2(n) 

+ (2880 + 11520 + 17280 + 11520 + 2880 n^) S-2,i,i{^)) (6 



Note that we reduced the number of different harmonic sums to 6: 

5_i(n), Si(n), 5_2,i(n), S2,-2(n), 52,i(n), S_2,i,i(n) 

by using in addition the differential operator D = ^ . Finally we can apply half-integer 
relation (see Chapter 4) with the function call ReduceToBasisfS", UseDifferentiation — )• 
True, UseHalflnteger True]: 

^NS,2,CaCfNf = 

— J (720 - 60 n - 1810 - 2085 - 8015 n"^ - 480 (-1)" - 9275 

135n4(l + n)^^ ^ ' 

- 1920 (-1)" + 1735 + 5400 n7 + 1350 - 180 ra^ C2 + 240 C2 - 1360 n'^ C2 
-10840 C2 - 19080 (2 - 13360 C2 - 3340 C2 + 1224 C2^ + 4896 C,2^ 
+7344 C2^ + 4896 C2^ + 1224 C2^ + 270 Cs + 2430 Cs + 7290 Cs 
+9990 Cs + 6480 Cs + 1620 Cs + (I8O + 1140 + 2940 + 3780 + 2400 
+600 n^) D [Half [Si (n),n] - 5i(n),n,2] + (240 n'' + 960 + 1440 + 960 
+240 n^) D [Half [Si (n),n] - Si(n),n,3] + (3340 n'' + 13360 + 20040 
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+13360 n'^ + 3340 - 720 n'^ (2 - 2880 (2 - 4320 (2 - 2880 n'^ (2 - 720 C2) 
B[Si{n),n,l] + (-360n^ - 1440n^ - 2160n6 - 1440n'^ - 360n^) B[Si{n),n,lf 
+D[Half[5'i(ra),n] - Si{n),n, 1] (-360 ra^ + 480 ra^ + 3960 + 5040 + 1920 + 1440 
+5760 C2 + 8640 C2 + 5760 (2 + 1440 C2 + (- 1440 - 5760 - 8640 
-5760 - 1440 n^) D[Si(n), n, 1]) + (270 + 2070 n'^ + 5850 + 7830 
+5040 + 1260 n^) D[Si(n), n, 2] + (360 n"^ + 1440 + 2160 + 1440 + 360 n^) 
B[Si{n),n,3] + (-720 - 2880 - 4320 - 2880 - 720 n^) B[S2,i{n),n, 1] 
+ (720 - 240 - 7440 - 13680 - 9600 - 2400 n^) 5_2,i(n) + Si{n) (180 n 
+540 + 720 + 3170 n'^ + 720 (-1)" n"^ + 9080 + 720 (-l)"" + 12540 
+8360 + 2090 - 1200 C2 - 4800 C2 - 7200 C2 - 4800 C2 - 1200 C2 
-540 n'' Cs - 2160 Cs - 3240 Cs - 2160 n'^ Cs - 540 Cs + (-2400 - 9600 
- 14400 -9600 - 2400n^) D[Half[S'i(n), n] - Si{n),n,l] + (-360^"^ - 1440n^ 
-2160n^ - 1440n^ - 360n®) D[Half[S'i(n), n] - 51(71), n, 2] + (-540n'' - 2160n^ 
-3240 - 2160 - 540 n*^) B[Si{n),n,2] + (-1440 71"^ - 5760 - 8640 - 5760 
-1440n^) S'_2,i(n)) + (720 n^ + 2880 n^ + 4320 n^ + 2880 n^ + 720 n^) S'2-2(n) 
+ (2880 n'' + 11520^5 + 17280 n*^ + 11520 + 2880 n^) 5_2,i,i(n)) 

In this way only the 5 multiple harmonic sums 

Si{n),S-2,i{^),S2-2{n),S2,i{n) and S'_2,i,i(n) 



together with D and the half-integer function Half remain. Hence we were able to 
reduce the number of contributing sums from 15 in (6.2) to 5 in (6.5). 
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